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A  problem  of  routing  a  U.S.  Coast  Guard  buoytender  to 
service  aids  to  navigation  is  formulated  as  a  symmetric 
traveling  salesman  problem.  A  heuristic  algoritnm  is  devel¬ 
oped  which  seeks  the  minimum  distance  tour  which  can  be 
taken  by  the  buoytender  to  visit  the  aids  to  navigation.  A 
user's  guide  is  provided. 

The  algorithm  is  programmed  in  Convergent  Technologies 
FORTRAN  for  use  on  the  Coast  Guard  Standard  Terminal. 
Several  problems  are  solved  by  the  algorithm  producing  solu¬ 
tions  that  are  optimal  or  nearly  optimal. 
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I.  IMTBODOCTION 


One  of  the  U.S.  Coast  Guard's  missions  is  maintaining 
aids  to  navigation  in  the  waters  of  the  United  States. 
These  aids  require  periodic  servicing  to  ensure  they  are  on 
station  and  they  are  showing  their  proper  watch  characteris¬ 
tics.  The  responsibility  of  maintaining  many  of  these  aids 
falls  on  the  Coast  Guard's  fleet  of  buoy tenders.  A  buoy- 
tender  in  its  area  of  responsibility  may  have  as  many  as  200 
aids  which  it  is  required  to  maintain. 

The  buoytender  has  scheduled  ATON  runs  (Aids  TO 
Navigation)  and  emergency  outages  it  must  handle  during  a 
fiscal  year.  An  emergency  outage  occurs  when  an  aid  is 
reported  by  a  mariner  to  be  showing  improper  characteris¬ 
tics  or  is  missing  from  its  station.  These  outages  are 
handled  shortly  after  they  are  reported.  Several  times 
during  the  year  the  buoytender  plans  an  ATON  run  to  handle 
the  periodic  servicing  of  aids,  to  ensure  that  they  maintain 
their  proper  watch  characteristics. 

Prior  to  a  scheduled  run  the  tender  will  raceive  SANDS 
forms  from  the  unit's  district  office  for  each  of  the  aids 
which  are  due  for  service  or  relief.  SANDS  is  a  database 
where  information  is  kept  on  a  district's  navigational  aids. 
The  district's  aids  to  navigation  branch  identifies  aids 
which  are  due  for  service  or  relief  and  notifies  the  tender. 
The  district  provides  a  list  of  aids  which  need  to  be 
serviced  during  the  scheduled  run  and  of  work  needed  on  each 
aid.  The  tender  then  makes  plans  accordingly  to  schedule 
and  complete  the  work. 

The  number  of  aids  serviced  on  a  particular  scheduled 
ATON  run  may  be  as  few  as  5  (a  day  of  local  work)  to  as  many 
as  70  (an  extended  four-week  trip  servicing  the  west  coast 
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of  Alaska  from  the  Aleutian  Islands  to  Point  Hope).  Ihe 
maintenance  of  the  floating  aids  and  fixed  structures 
entails  the  travel  of  the  tender  from  aid  to  aid.  The  tend¬ 
er’s  Operations  Officer  usually  plans  the  initial  route  to 
he  travelled  to  service  the  aids  due  for  routine  maintenance 
and  presents  his  proposal  to  the  Commanding  Officer,  who  may 
then  modify  this  proposal. 

There  are  several  factors  which  go  into  the  decision 
process  to  produce  the  desired  route.  Two  of  these  factors 
are  (1)  which  aids  are  scheduled  for  service,  and  (2)  the 
distances  between  the  aids.  The  desired  route  is  usually 
the  route  of  shortest  distance  which  visits  all  the  sched¬ 
uled  aids  and  returns  the  tender  to  her  point  of  origin. 
This  problem  is  a  classic  problem  which  operations  analysts 
normally  call  'the  traveling  salesman  problem.'  Present 
route  selection  methods  are  (1)  traditional  or  previously 
followed  routes  and  (2)  routes  created  by  Operations 
Department  personnel  sitting  down  with  chart  aud  rule  and 
selecting  a  route. 

It  is  the  intent  of  this  thesis  to  propose  and  implement 
a  computer  assisted  approach  as  an  alternative  method  for 
the  tuoytender  route  selection  problem.  The  succeeding 
chapters  of  this  paper  contain  a  discussion  of  the  trav¬ 
eling  salesman  problem  and  how  it  relates  to  the  buoytender 
problem,  a  brief  survey  of  possible  solution  methods  to  the 
problem,  a  description  of  the  solution  process  of  the 
selected  method,  and  the  results  of  the  route  selection 
program  with  conclusions. 
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II-  NATURE  OF  JHE  £ROBLEM 


In  this  chapter  the  principles  of  the  traveling  salesman 
problem  will  be  presented  and  it  will  be  shown  how  the  buoy- 
tender  problem  can  be  described  as  a  traveling  salesman 
problem.  The  discussion  will  include  a  general  description, 
definitions  about  solution  methods,  and  the  solution 
requirements  of  the  problem.  Further,  a  description  of  the 
computing  hardware  on  which  the  selected  solution  method  is 
programmed  is  provided.  This  provides  an  idea  of  the  capa¬ 
bilities  of  the  hardware  since  the  hardware  has  a  bearing  on 
the  selection  of  the  solution  method. 

A.  DESCRIPTION  OF  Tfi£  TRAVELING  SALESMAN  PRO BLEB 

A  brief  description  of  the  traveling  salesman  problem 
(or  TSP)  is  as  follows.  A  salesman  has  n  cities  he  must 
visit.  The  salesman  starts  in  one  of  the  cities  and  must 
travel  so  that  he  enters  and  departs  each  of  the  other 
cities  only  once.  Upon  visiting  the  last  city  the  salesman 
will  return  to  the  city  from  which  he  originally  started  his 
travels.  The  desired  solution  to  the  problem  is  the  route 
the  salesman  can  traverse  which  is  of  the  shortest  total 
distance  possible. 

It  is  usually  assumed  that  a  TSP  can  be  associated  with 
some  n  by  n  distance  matrix.  The  elements  of  the  matrix, 
d(*  ,  are  distances  from  city  i  to  city  j,  (i  =  1,...,n), 
(j=1, ...  ,n)  ,and  where  dji=oo,  i=1,...,n.  The  traveling 
salesman  is  not  allowed  to  leave  city  i  and  to  return  to 
city  i  in  his  tour. 

A  T5E  can  be  associated  with  a  graph,  also.  A  city 
which  the  salesman  must  visit  is  represented  by  a  node  in 


this  graph.  There  is  an  arc  (i,j)  in  the  graph  with  length 
dt'j  for  each  node. 

A  distance  matrix  may  be  classified  as  sparse  or  dense. 
A  sparse  matrix  will  have  a  majority  of  its  entries  egual  to 
infinity.  The  corresponding  graph  is  sparse  since  it  has 
few  edges.  A  dense  matrix  is  a  matrix  where  a  majority  if 
not  all  its  entries,  with  the  exception  of  the  d pp  *s,  are 
less  than  infinity.  The  number  of  possible  permutations  of 
cities  for  a  dense  matrix  is  of  order  (n— 1)1.  Computation 
time  and  memory  required  for  solution  by  total  enumeration 
grows  rapidly  with  the  size  of  n  due  to  the  large  number  of 
possible  solutions,  each  permutation  of  cities  being  a 
feasible  tour  and  a  possible  solution. 

The  distance  matrix  used  by  the  buoytender  for  route 
selection  will  be  a  dense  matrix.  With  the  large  number  of 
permutations  or  possible  solutions  to  the  problem,  some 
means  other  than  total  enumeration  must  be  used  to  solve  the 
buoytender  problem. 

The  traveling  salesman  problem  may  be  broken  into  two 
classes,  symmetric  and  asymmetric.  In  a  symmetric  TSP  the 
d;j  =d a  for  all  i^j,  while  in  an  asymmetric  TSP  this  need 
not  be  true.  Further,  a  TSP  may  or  may  not  be  required  to 
satisfy  the  triangle  inequality.  Consider  three  cities  i, 
j,  k.  The  triangle  inequality  states  that 


dpi  <  dik  ♦  d»<: 


for  all  i^j^k^i. 


The  conditions  which  make  a  TSP  asymmetric,  the  unequal 
distances  between  cities,  may  cause  the  problem  to  violate 
the  triangle  inequality. 

The  TSP,  both  symmetric  and  asymmetric,  belongs  to  the 
N?  Complete  class  of  problems  [Bef.  1,2].  This  has  one 
major  implication  with  respect  to  the  complexity  of  solving 
the  TSP:  no  polynomial-time  algorithms  are  known  or  seem 
likely  tc  be  devised  for  exact  solution  of  the  TSP.  Only 
exponential-time  exact  algorithms  are  known. 


There  have  been  many  algorithms  proposed  for  the  solu¬ 
tion  of  the  TSP,  each  algorithm  showing  some  advantages  and 
disadvantages  to  its  method  of  attacking  and  solving  the 
problem.  Several  of  the  published  algorithms  are  tailored 
to  a  particular  size  of  TSP.  The  TSP  may  be  categorized  by 
size  as  either  small  or  large.  A  reference  could  not  be 
found  in  the  literature  which  would  give  a  hard  numerical 
value  by  which  the  size  of  a  TSP  could  be  judged,  but  a 
small  TSE  is  generally  considered  to  be  less  then  15  cities. 
A  large  ISP  is  often  considered  to  be  in  the  neighborhood  of 
45  or  more  cities.  There  are  solution  methods  which  solve 
the  TSP  exactly  for  small  up  to  large  problems.  There  are 
some  methods  which  can  solve  some  large  problems  exactly, 
but  as  a  rule  large  problems  are  solved  with  heuristics  or 
approximate  methods. 

The  circumstances  of  the  buoytender  problem  satisfy  the 
conditions  of  a  class  of  the  traveling  salesman  problem. 
The  buoytender  problem  can  be  characterized  as  a  symmetric 
TSP,  i.e.,  the  distance  traveled  between  any  two  aids  is  the 
same  regardless  cf  direction  travelled  by  the  tender.  The 
size  cf  buoytender  problem  will  vary  with  the  number  of  aids 
scheduled  to  be  serviced  on  a  specific  ATON  run;  therefore 
the  solution  method  for  the  buoytender  problem  will  have  to 
handle  both  small  and  large  problems. 

The  triangle  inequality  as  applied  to  the  buoytender 
problem  holds  true,  tut  how  it  holds  needs  some  explanation, 
with  all  distances  positive,  the  triangle  inequality  states 
that  given  three  points  A,  B,  and  C,  the  distance  from  A  to 
C  must  be  less  than  or  equal  to  the  distance  from  A  to  B 
plus  the  distance  from  B  to  C  (see  Figure  2.1)  .  In  the 
buoytender* s  problem  the  straight-line  distance  from  A  to  C 
may  be  less  than  (A  to  3)  ♦  (B  to  C),  but  some  situations 
will  require  that  A  to  C  equal  (A  to  B)  ♦  (B  to  C)  (see 
Figure  2.1), 


d4C  =  dA6  *  dBC  • 


)  '  6 
J  / 


UhIc^ 


Figure  2.1  Triangle  Ineguality  /  Obstructed  A  to  C  Path 

The  buoytender  may  be  unable  to  take  the  direct  route  from  A 
to  C  due  to  shoal  water  or  an  obstruction.  To  get  from  A  to 
C  the  tender  will  have  to  take  the  routes  A  to  B  and  B  to  C. 
This  is  not  an  unusual  situation  for  buoytenders  since  the 
aids  being  serviced  often  mark  shoal  water  or  obstructions. 


B.  DESCRIPTION  OF  TEE  SOLUTION  REQUIREMENTS 

There  are  many  proposed  solution  methods  to  the  TSP  and 
they  may  be  divided  into  two  groups:  exact  methods  and 
approximate  or  heuristic  methods.  An  exact  method  or  algo¬ 
rithm  has  the  property  that  upon  termination  of  the  algo¬ 
rithm,  the  user  will  have  the  best  possible  or  global 
optimum  solution.  Heuristics  or  approximate  methods,  on  the 
other  hand,  terminate  with  a  local  optimum  solution  or  the 
best  solution  found  thus  far.  This  local  optimum  solution 
may  be  the  global  optimum  but  usually  there  is  no  proof  that 
the  global  optimum  has  been  found.  Heuristics  or  approxi¬ 
mate  methods  solve  by  checking  a  subset  of  the  numerous 
possible  solutions  and  terminate  when  the  subset  has  teen 
searched  for  the  best  answer  it  contains,  a  specified  time 


limit  for  the  run  has  beer,  reached,  or  the  memory  capacity 
of  the  computing  machine  has  been  reached. 

A  further  consideration  in  selecting  a  method  to  solve 
the  TSP  is  that  there  should  be  very  little  interactive  work 
for  the  user.  The  user  interactive  work  necessary  to 
operate  the  program  should  be  understandable  to  a  majority 
of  the  intended  users,  who  in  this  case  are  the  buoytender’s 
Operations  Officer  or  Navigator,  and  Quartermaster 
personnel.  By  requiring  the  interactive  work  to  be  minimal 
and  simple  in  nature  the  possiblity  of  operator  introduced 
error  will  be  reduced.  Further,  if  the  interactive  work  is 
of  a  simple  nature,  the  program  will  be  useable  by  a 
majority  of  the  intended  users  without  extensive  training. 

Although  it  is  desired  to  obtain  the  global  optimum,  it 
is  not  mandatory.  The  computer  route  selection  method  is  to 
be  a  tool  used  in  planning,  and  the  exact  solution  is  not 
critical  to  the  operation  of  the  buoytender.  By  using 
heuristics,  computer  time  and  memory  may  be  reduced  since  a 
heuristic  generally  searches  only  a  subset  of  all  possible 
solutions.  If  a  computer  solution  can  be  produced  which  is 
tetter  than  the  traditional  route  or  a  route  produced  by 
hand  calculation,  then  a  goal  has  been  met,  that  goal  being 
reducing  the  total  distance  traveled  by  the  tender  thereby 
saving  fuel  and  time.  The  computer  solution  is  then  a 
viable  alternative  tc  be  added  to  the  present  methods  of 
route  selection. 

C.  CCHPOTING  HAHDHAEI 


The  computer  on  which  the  selected  solution  method  is 
programmed  is  the  Coast  Guard’s  C3  Standard  Terminal.  The 
reason  for  selecting  this  computer  is  that  it  is  presently 
being  installed  throughout  the  Coast  Guard  as  a  primary 
unit-level  computer.  Most  intend  users  should  have  access 
to  a  C3  Terminal. 


The  C3  Standard  Terminal  is  a  standalone  16  bit  micro¬ 
computer.  The  terminal  comes  in  two  basic  configurations, 
the  Integrated  Work  Station  (IWS)  and  the  Application  Work 
Station  (AWS)  .  The  IWS  is  a  master  or  standalone  unit 
equipped  with  a  hard  disk  and  8  inch  floppy  mass  storage 
system,  a  300/1200  baud  modem,  a  printer,  and  one  of  three 
RAM  memory  configurations  -  256k,  384k,  or  512k.  The  IWS 
may  be  programmmed  in  BASIC,  FORTRAN,  COBOL,  Pascal,  and 
Assembly.  The  AWS  is  more  limited  in  its  capabilities  than 
the  IWS.  The  AWS  is  limited  to  256k  RAM  and  must  be 
networked  to  a  IWS  station  to  access  the  peripherals.  With 
a  program  in  RAM  the  AWS  will  operate  as  a  standalone 
computer,  but  does  not  have  the  power  or  computational  speed 
of  the  IWS. 

The  solution  method  selected  for  programming  will  have 
to  perform  within  the  limitations  of  the  C3's  memory.  The 
C3's  RAM  memory  will  be  a  critical  factor  in  selection  of 
the  solution  method.  Since  the  computer  is  a  purchased 
standalone  and  not  a  leased  system,  CPU  time  should  not  be  a 
critical  factor.  It  is  still  desired,  though,  to  keep  CPU 
time  from  becoming  excessive  since  the  computer  is  also 
needed  for  ether  work.  The  computer  could  be  allowed  to  run 
overnight  (off  duty  hours)  for  large  problems,  and  therefore 
a  solution  method  requiring  more  than  12  hours  to  provide  an 
acceptable  solution  would  be  excessive.  A  solution  method 
which  provides  an  answer  in  under  an  hour  could  be  run 
during  duty  hours  with  very  little  impact  on  other  users, 
particularly  if  it  was  run  on  a  slave  terminal. 

In  the  next  chapter  we  take  a  brief  look  at  some 
suggested  solution  methods  for  the  TSP  and  their  applica¬ 
bility  to  the  buoytender  problem. 


III.  A  Mlil  SOHVEX  OF  SOLUTION  HSTHOpS  FOE  IHJ  ISP 


In  this  chapter  a  brief  survey  of  the  general  solution 
methods  for  the  TSP  and  some  of  their  associated  algorithms 
will  be  presented.  The  number  of  algorithms  proposed  for 
solving  the  TSP  is  extensive.  To  present  all  of  the  avail¬ 
able  algorithms  or  to  discuss  them  in  detail  is  worthy  of  a 
dissertation  in  itself.  This  survey  of  solution  methods  is 
made  because  a  combination  of  a  couple  of  the  methods  will 
be  used  to  solve  the  buoytender  problem.  This  survey  will 
be  brief  in  nature  presenting  the  general  solution  methods 
and  briefly  describing  some  of  the  more  popular  of  published 
algorithms  which  fall  under  these  solution  method  headings.1 

The  solution  methods  may  be  classified  under  four 
general  headings.  These  headings  are  Tour  Building,  Subtour 
Elimination,  The  One-Tree  Formulation,  and  Tour  Improvement. 
These  headings  describe  how  the  traveling  salesman  problem 
is  approached. 

A.  TCDH  BUILDING 

Tcur  Building  solution  methods  to  construct  a  tour  using 
the  n  nodes  and  the  available  arcs  in  the  problem.  Tour 
Building  contains  both  exact  and  approximate  algorithms. 
Some  of  the  approaches  taken  by  Tour  Building  Algorithms  to 
solve  the  TSP  are  dynamic  programming,  'branch  and  bound* 
methods,  and  tour  construction  heuristics. 


iThree  excellent  articles  [Ref.  1#3,4]  have  been 
published  which  describe  the  TSP,  discuss  general  solution 
methods,  and  review  many  of  the  algorithms  which  solve  it. 
These  three  articles  review  the  research  conducted  on  the 
TSP  from  the  mid  1950’s  up  to  1983. 
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Held  and  Karp  [Bef.  5]  present  a  dynamic  programming 
approach  which  is  suitable  for  small  problems  and  can  be 
modified  to  an  approximate  method  for  larger  proDlems.  The 
dynamic  programming  approach  rapidly  consumes  available 
memory  as  the  size  of  the  problem  increases.  Due  to  memory 
limitations,  the  algorithm  Held  and  Karp  presented  is  exact 
only  up  to  13  cities.  They  also  present  a  dynamic  program¬ 
ming  method  which  approximates  the  solution  for  larger  prob¬ 
lems.  The  approximation  is  done  by  partitioning  the  large 
problem  into  several  smaller  problems  or  subproblems,  each 
subproblem  containing  13  or  fewer  nodes.  Each  partitioned 
subproblem  is  solved  for  its  optimal  tour  resulting  in  a  set 
of  subtours  for  the  original  problem.  The  subtours  are  then 
broken  and  linked  to  each  other  to  provide  a  tour  for  the 
original  large  problem. 

Little  et  al.  [Bef.  6]  present  a  branch  and  bound 
method.  The  algorithm  branches  on  whether  a  particular  arc 
is  included  or  excluded  from  the  tour  being  constructed. 
This  branching  partitions  the  solution  set  containing  all 
tours  into  subsets.  Each  arc  in  the  total  graph  can  be 
associated  with  s^me  subset  of  tours,  and  the  idea  is  to 
find  the  subset  containing  the  optimal  tour.  For  each 
subset  of  tours,  a  lower  bound  is  computed  for  the  tours 
within  the  subset.  The  tour  is  constructed  as  arcs  with 
favorable  lcwer  bounds  are  included  in  the  tour.  The  algo¬ 
rithm  continues  branching  and  computing  lower  bojnds  until  a 
subset  is  found  which  contains  one  tour  and  this  tour's 
distance  is  less  than  or  equal  to  the  lower  bounds  of  the 
other  subsets  of  tours.  This  tour  will  be  the  optimal  solu¬ 
tion  to  the  TSP.  The  memory  required  for  this  technique  may 
be  extensive  and  computing  time  grows  exponentially  with  the 
size  cf  the  problem. 

Many  tour  construction  heuristics  have  been  proposed. 
Bosenkrantz,  Stearns,  and  Lewis  [Bef.  2]  and  Golden  et  al. 


[Ref.  4]  review  several  of  these  heuristics.  The  two  most 
popular  forms  of  these  heuristics  are  the  nearest  neighbor 
algorithm  and  the  insertion  techniques. 

The  nearest  neighbor  algorithm  operates  exactly  as  its 
name  implies.  A  starting  node  is  chosen  and  its  nearest 
neighboring  node  is  selected  as  the  next  node  in  the  tour 
being  constructed.  The  algorithm  then  iteratively  selects  a 
node  not  yet  in  the  tour  which  is  nearest  to  the  previously 
selected  node  to  enter  the  tour.  This  selection  continues 
until  all  nodes  have  been  selected  and  then  the  last  node 
selected  is  connected  to  the  first  node  to  complete  the 
tour. 

The  insertion  techniques  include  the  nearest,  farthest, 
cheapest,  and  arbitrary  insertion  algorithms.  Insertion 
techniques  begin  with  a  starting  node  and  select  the  nearest 
node  to  create  a  sub  tour.  The  algorithm  iterates  selecting 
a  node  which  is  nearest  to  or  farthest  from  any  other  node 
in  the  subtour,  or  arbitrarily  selects  the  next  node  to 
enter  the  tour.  This  selected  node  is  then  inserted  into 
the  subtcur  wherever  it  provides  the  least  increase  (or  in 
the  case  of  the  farthest,  the  greatest  increase)  in  distance 
in  the  new  subtour.  These  heuristics  generally  provide  a 
suboptimal  solution  to  the  TSP,  but  they  have  the  advantage 
of  being  very  simple  to  implement  and  have  very  fast  solu¬ 
tion  times  with  linimal  memory  consumption. 

B.  SOBTCOR  ELIHIHATICH 

The  Subtour  Elimination  solution  method  [Ref.  t]  attacks 
the  TSP  by  solving  an  n  by  n  assignment  problem  with  the 
added  constraints  that  the  final  solution  must  contain  a 
cycle  and  the  cycle  cannot  be  of  size  n-1  or  less.  The 
number  cf  integer  constraints  necessary  to  eliminate  all 
subtours  or  cycles  of  size  n-1  or  less  is  enormous. 
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Therefore,  the  initial  solution  is  generally  found  using  a 
'relaxed*  form  of  the  assignment  problem.  The  relaxed 
assignment  problem  emits  the  subtour  constraints  in  its 
formulation. 

Mathematically  stated  the  relaxed  assignment  problem  is 


Min 
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xq  =0 

or  x2i 

=  1 

The  solution  to  the  relaxed  assignment  problem  provides 
the  initial  lover  bound  on  the  optimal  tour  length.  If  the 
optimal  assignment  solution  is  a  tour  then  it  is  an  optimal 
solution  to  the  TSP.  If  the  assignment  solution  is  not  a 
tour  then  there  exist  subtours  which  must  be  removed  until  a 
single  tour  exists.  If  a  subtour  exists  with  k  arcs  in  the 
subtour,  then  there  are  k  possible  subproblems  to  which  the 
problem  may  branch.  These  k  subproblems  each  have  an  addi¬ 
tional  constraint,  the  constraint  for  subprofclem  i 
(i=1,...,k)  being  the  exclusion  of  arc  i  from  the  problem  to 
eliminate  the  subtour.  Instead  of  branching  into  k  subprob¬ 
lems  the  approach  is  to  branch  into  two  subproblems.  An  arc 
is  selected  from  the  subtour  and  the  exclusion  of  this  arc 
becomes  a  new  constraint  in  one  of  the  subproblems  while 
inclusion  of  the  arc  in  the  tour  becomes  a  new  constraint 
for  the  alternate  subproblem.  The  modified  assignment 
problem  is  solved  again  for  each  of  these  new  constraints 
and  provides  the  lower  bound  for  its  respective  branch.  In 
this  way  only  the  necessary  subtour  elimination  constraints 
are  added  as  needed  to  the  assignment  problem,  rather  than 
attempting  to  add  all  possible  subtour  constraints  to  the 
initial  assignment  problem.  Ideally  the  branch  solution 
subset  with  the  lowest  bound  is  then  selected  to  he  tested 
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to  see  if  it  is  a  tour,  and,  if  it  is  not,  the  problem 
tranches  again  on  the  lowest  bound.  This  branching 
continues  until  all  subtours  have  been  eliminated  and  a  tour 
is  found  which  has  a  distance  less  than  or  equal  to  the 
lower  bounds  of  the  ether  solution  subsets. 

There  are  two  basic  methods  for  searching  the  tree 
created  by  the  branch  and  bound  process.  One  method  is  a 
breadth-first  search.  As  the  program  is  branching  and 
computing  bounds,  a  tree  is  created  with  leaves  or  terminal 
nodes.  A  breadth- first  search  requires  storage  for  the 
solutiens  of  each  leaf  and  a  search  through  these  leaves  for 
the  next  subproblem  on  which  to  branch.  This  requires 
extensive,  usually  exponential,  storage.  Branching  at  the 
lowest  bound  requires  a  breadth-first  search.  The  second 
method  for  searching  the  tree  is  a  depth-first  search.  The 
depth-first  search  is  a  more  localized  approach  to  handling 
♦he  search.  The  depth-first  search  need  only  store  the 
immediate  solution  while  searching  a  branch.  If  during  the 
search  it  is  found  that  no  further  branching  is  possible 
from  the  current  node,  the  process  must  return  to  a  previous 
leaf  to  continue  the  search  for  the  optimal  tour.  The 
depth-first  search  backtracks  and  recreates  the  previous 
solution  at  that  leaf.  The  depth-first  search  does  not 
require  exponential  storage  like  the  breadth-first  search, 
but  it  usually  requires  more  computation  time. 

The  subtour  elimination  method  is  described  as  being 
exact.  The  subtour  elimination  method  usually  requires  less 
computation  time  than  the  branch  and  bound  approach  of 
little  et  al.  since  the  bounds  obtained  with  the  assignment 
relaxation  are  tighter. 


C.  TEE  CNE-TBEE  FORHULATIOH 


The  one-tree  formulation  was  introduced  and  developed  by 
Held  and  Karp  [Ref.  7,8].  The  general  method  for  the  solu¬ 
tion  of  the  TSP  begins  with  the  construction  of  a  minimum- 
weight  ' cne-tree.'  If  the  one-tree  is  not  a  tour  then  an 
integer  linear  program  may  be  used  to  obtain  the  final  solu¬ 
tion. 

To  assist  in  the  description  of  this  solution  method 
several  terms  need  defining.  First  a  tree  is  defined  as  a 
connected  graph  without  cycles.  A  minimum- weight  spanning 
tree  is  a  tree  with  minimum  total  weight  on  its  edges.  A 
one- tree  is  a  tree  to  which  one  edge  has  been  added  yielding 
exactly  one  cycle.  A  minimum-weight  one- tree  is  a  one-tree 
with  minimum  total  weight  on  its  edges.  Finding  a  minimum- 
weight  one-tree  is  a  relaxed  version  of  a  TSP  since  the 
solution  to  a  TSP  is  a  min imum-weight  one-tree  having  every 
vertex  of  degree  two. 

In  the  solution  method  proposed  by  Held  and  Karp 
[Ref.  7]  the  algorithm  begins  by  constructing  a  minimum- 
weight  one- tree.  A  minimum-weight  one-tree  can  be  can  be 
found  by  first  constructing  a  minimum- weight  spanning  tree 
through  nodes  2  to  n  and  then  adding  to  the  graph  the  two 
arcs  of  least  weight  from  node  1  [Ref.  9],  A  minimum- 
weight  one-tree  may  also  be  constructed  by  constructing  a 
minimum- weight  spanning  tree  through  nodes  1  to  n  and  then 
add  to  the  graph  the  next  minimum-weight  edge  not  yet  used 
from  the  distance  matrix.  If  the  one-tree  constructed  is  a 
tour  (i.e.  the  one-tree’s  vertices  are  all  of  degree  two) 
then  the  TSP  is  solved.  Otherwise,  the  one-tree  must  be 
converted  into  a  tour. 

Integer  linear  programming  is  used  to  obtain  an  optimal 
tour  if  the  one-tree  formulation  does  not  satisfy  a  tour. 
Held  and  Karp  introduce  the  concept  of  a  ’gap'  function  and 


then  use  a  special  integer  linear  program  to  minimize  this 
function.  In  essence  there  is  a  transformation  on  the  one- 
tree  variables  which  allows  the  integer  linear  program  to 
search  for  the  optimal  tour.  Held  and  Karp  suggested  two 
methods  for  optimizing  the  integer  linear  program,  an  ascent 
method  and  a  branch  and  bound  method  with  the  ascent  method 
embedded  in  it.  Other  improvements  to  the  Held  and  Karp 
method  were  suggested  by  Hansen,  Krarup  [Ref.  10]  and  Houck, 
Picard,  Vemuganti  [Ref.  11].  Bazarra  and  Goode  [Bef.  12] 
did  further  work  on  Held  and  Karp's  approach  using  optimiza¬ 
tion  of  a  lagrangian  dual  in  lieu  of  solving  integer  linear 
program.  They  proposed  a  branch  and  bound  scheme  with  subg¬ 
radient  optimization  of  the  dual  to  transform  the  one-tree 
into  the  optimal  tour.  The  one-tree  method  as  described  is 
exact  and  may  require  extensive  computation  time. 

0.  TCOR  IMPROVEMENT 

The  Tour  Improvement  method  assumes  that  an  arbitrary 
tour  is  available.  The  method  operates  by  perturbing  or 
exchanging  arcs  in  the  tour  until  a  better  tour  is  found. 
The  Tour  Improvement  method  terminates  when  a  better  tour 
cannot  be  found.  In  its  simplest  form,  it  is  possible  for 
this  method  to  check  (n-1)  1/2  arc  exchanges  for  the  symme¬ 
tric  TSP  before  terminating  (i.e.,  it  investigates  all 
available  answers)  . 

Dantzig,  Fulkerson,  and  Johnson  [Ref.  13,14]  proposed  a 
algorithm  which  starts  with  an  arbitrary  tour  and  uses 
integer  linear  programming  to  improve  the  tour.  The  TSP  is 
transformed  into  an  integer  linear  program  and  solved  using 
the  simplex  method.  Dantzig,  Fulkerson,  and  Johnson 
[Bef.  15]  mention  that  the  number  of  constraints  needed  to 
characterize  the  problem  is  'astronomical. '  Instead  of  using 
all  the  constraints,  they  begin  by  using  a  relaxed  version 


of  the  problem  and  -dd  constraints  as  needed  to  the  integer 
linear  program  to  maintain  feasibility.  In  their  method 
infeasibility  appears  as  'loop*  or  subtour.  To  continue 
solving  the  problem  a  constraint  is  added  which  removes  the 
subtcur  yet  does  not  eliminate  any  of  the  available  tours. 
This  algorithm  requires  the  addition  of  constraints  during 
the  solution  process  to  maintain  feasibility.  Several 
researchers  have  suggested  improvements  on  Dantzig, 
Fulkerson,  and  Johnson's  method  for  solution  of  the  TSP. 
[Ref.  16,17,18,19] 

Croes  [Ref.  20]  proposed  a  heuristic  whereby  a  subset  of 
all  possible  transformations  is  tested,  a  transformation 
being  the  transformation  of  one  tour  into  another  tour.  He 
called  these  transformations  'inversions'  because  they 
inverted  the  sequence  of  nodes  in  part  of  the  tour  to  create 
a  possible  improvement.  Lin  [Ref.  21]  later  describes  Croes 
'inversion  free  tours'  as  2-optimal  tours  and  goes  on  to 
describe  k-optimality  where  k  is  some  fixed  number  less  than 
n.  The  general  idea  is  to  exchange  k  arcs  iteratively  in 
the  tour  while  testing  for  improvements.  Lin  and  Kernighan 
[Ref.  22]  presents  a  modified  k-opt  algorithm  where  k  is  not 
fixed.  In  this  algorithm  k  may  vary  from  2  to  n.  The 
floating  k-opt  algorithm  is  also  described  by  Christofides 
and  Eilon.  [Ref.  23] 

The  advantage  of  the  'inversions'  or  k-opt  solution 
method  is  that  for  a  given  problem,  the  memory  needed  to 
solve  the  problem  is  fixed.  There  are  no  constraints  which 
must  be  added  to  the  problem,  and  the  decision  rule  used  is 
simple.  The  k-opt  method  is  exact,  but  may  be  terminated 
before  optimality  is  reached  to  provide  a  satisfactory 
(suboptimal)  solution.  This  generally  reduces  run  time  and 
still  provides  a  gocd  solution  to  the  problem.  The  k-opt 
method  may  also  be  programmed  to  require  very  little  user 
interaction;  in  the  best  case  the  user  need  only  enter  the 
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distance  matrix  for  the  nodes  desired  in  the  tour.  The 
k-opt  algorithm  may  require  extensive  computation  time, 
particularly  if  during  solution  k  is  found  to  he  greater 
than  2  on  many  of  the  iterations. 

This  chapter  has  presented  a  general  survey  of  some  the 
solution  methods  for  the  TSP  and  some  algorithms  associated 
with  these  solution  methods.  The  next  chapter  will  discuss 
which  of  these  solution  methods  were  employed  to  solve  the 
bouytender  problem.  The  buoytender  routing  algorithm  will 
then  be  presented  and  discussed  in  detail. 


IV.  BDOYTENDER  ALGORITHM 

This  chapter  will  discuss  the  algorithm  with  which  we 
propose  to  solve  the  buoytender  problem.  The  reasons  for 
the  algorithm’s  selection  and  some  background  on  its  devel¬ 
opment  will  be  presented,  together  with  a  detailed  descrip¬ 
tion  of  how  the  algorithm  addresses  the  buoytender  problem. 

A.  A1G0BITHM  SELECTION 

The  algorithm  is  a  combination  of  two  of  the  previously 
mentioned  solution  methods  in  that  it  combines  a  Tcur 


Building  heuristic  with  a  Tour  to  Tour  Improvement 
heuristic.  The  concept  of  combining  a  Tour  Building 
heuristic  with  a  2-cpt  or  3-opt  heuristic  was  proposed  by 
Golden  et  al.  [Bef.  24]  as  a  relatively  fast  computational 
solution  method  which  would  provide  an  optimal  or  near 
optimal  tour.  Their  'composite*  algorithm  is  the  foundation 
on  which  the  buoytender  algorithm  was  developed. 

The  heuristics  used  in  the  buoytender  algorithm  are  the 
nearest  neighbor  algorithm  and  the  k-opt  algorithm.  The 
nearest  neighbor  algorithm  is  used  to  construct  an  initial 
tour  while  a  version  of  the  k-opt  proposed  by  Lin  and 
Kernighan  £Ref.  22]  is  used  to  improve  this  initial  tour. 
The  nearest  neighbor  and  k-opt  heuristics  were  selected  for 
use  in  the  buoytender  algorithm  because  very  little  interac¬ 
tive  work  is  required,  the  interactive  work  is  simple,  no 
constraints  need  to  be  added  to  maintain  feasibility,  compu¬ 
tational  time  is  reasonable,  and  the  algorithm  operates  with 
a  fixed  amount  of  memory. 

The  buoytender  algorithm  as  presented  requires  very 
little  user  interaction.  The  nearest  neighbor  algorithm 


needs  only  the  entry  of  the  distance  matrix  for  the  problem 
and  the  selection  of  a  starting  node.  The  k-opt  algorithm 
needs  the  same  distance  matrix  and  an  initial  tour.  The 
heuristics*  decision  processes  will  provide  a  feasible  and 
possibly  optimal  solution  in  a  reasonable  period  of  time. 
The  nearest  neighbor  algorithm  was  selected  to  construct  the 
initial  tour  because  it  is  computationally  faster  than 
insertion  techniques,  and  provides  a  good  initial  solution 
to  the  problem.  [Ref.  25] 

A  critical  factor  in  solving  the  TSP  on  a  microcomputer 
is  the  memory  required  by  the  algorithm  used.  The  amount  of 
memory  required  for  solution  of  a  TSP  of  size  n  is  fixed  for 
the  nearest  neighbor  and  k-opt  heuristics.  With  known 
memory  requirements  for  the  heuristics,  a  program  can  be 
developed  to  fit  the  Coast  Guard  Standard  Terminal.  The 
nearest  neighbor  algorithm  requires  an  n  by  n  distance 
matrix  and  an  n  arrary  in  which  to  store  the  tour  as  it  is 
constructed.  The  k-cpt  algorithm  requires  the  same  storage 
as  the  nearest  neighbor  plus  a  n  by  n  decision  matrix,  and  a 
2n  array  for  recording  nodes  selected  for  inversion  in  the 
tour. 

B.  PROGRAMMING  THE  ALGORITHM 

The  tuoytender  algorithm  was  developed  as  a  series  of 
modules  which  were  then  programmed  as  subroutines  in  the 
operating  program.  The  algorithm  was  programmed  in  FORTRAN 
for  operation  on  the  Coast  Guard  Standard  Terminal.  FORTRAN 
was  selected  because  it  is  efficient  and  a  majority  of 
programmers  are  familiar  with  it. 

There  are  some  differences  between  the  heuristics  in  the 
literature  and  the  buoytender  algorithm.  The  nearest 
neighbor  algorithm  as  used  by  the  buoytender  algorithm  is 
the  same  as  outlined  by  Rosenkrantz,  Stearns,  and  lewis 


[Bef.  2],  The  k-opt  heuristic  used  in  the  buoytender  algo¬ 
rithm  is  a  modified  version  of  the  heuristic  presented  by 
Lin  and  Kernighan  £  Bef.  22].  There  are  three  major  Infer¬ 
ences  between  the  k-cpt  in  the  bouytender  algorithm  and  the 
k-opt  suggested  by  Lin  and  Kernighan.  First,  the  Lin  and 
Kernighan  k-opt  heuristic  selects  all  possible  nodes  which 
will  result  in  a  shorter  tour  before  creating  the  improved 
tour.  The  buoytender  algorithm  improves  the  tour  by 
exchanging  nodes  whenever  a  favorable  selection  is  found. 
This  improvement-as-ycu-go  procedure  results  in  a  simplifed 
selection  process  for  the  next  set  of  nodes  to  be  tested. 
Second,  the  Lin  and  Kernighan  k-opt  heuristic  has  an  added 
facility  for  'limited  backtracking'  used  when  a  particular 
exchange  gives  an  improvement  of  zero.  The  backtracking 
procedure  searches  until  a  gain  greater  than  zero  can  be 
found  which  improves  the  tour.  In  the  interest  of  short¬ 
ening  computational  time,  the  buoytender  algorithm  does  not 
'backtrack*  when  a  gain  of  zero  is  found,  but  instead  treats 
zero  gain  as  no  improvement  and  continues  with  the  next 
selection.  The  third  major  difference  between  the  Lin  and 
Kernighan  k-opt  heuristic  and  the  buoytender  algorithm  is 
related  to  what  Lin  and  Kernighan  call  'reduction.'  Their 
heuristic,  after  producing  several  locally  optimal  tours, 
checks  fcr  arcs  which  appear  in  each  of  the  tours.  These 
'good'  arcs  are  not  allowed  to  be  broken  in  further  computa¬ 
tions  for  other  locally  optimal  tours  thereby  reducing  the 
number  of  links  to  be  checked  for  improving  the  tour.  While 
decreasing  run  time,  this  procedure  requires  more  memory  and 
is  therefore  omitted  from  the  buoytender  algorithm. 

The  buoytender  algorithm  can  be  divided  into  three  basic 
operations.  First,  the  algorithm  selects  a  random  sample  of 
k  nodes,  k<n.  This  random  sample  of  nodes  becomes  the  set 
of  starting  nodes  for  the  next  major  operation  in  the  algo¬ 
rithm,  to  create  k  nearest  neighbor  tours.  The  third 


operation  is  to  use  the  k  nearest  neighbor  tours  as  initial 
tours  for  the  k-opt  improvement.  The  best  of  the  k-opt 
improvement  tours  is  presented  as  the  solution  to  the 
problem. 

A  further  modification  was  made  to  the  algorithm  based 
on  some  computational  results  obtained  from  use  of  an  early 
version  of  the  algorithm.  During  programming  of  the  buoy- 
tender  algorithm  several  benchmark  statistics  were  obtained 
from  the  program.  Two  of  these  statistics  were  the  initial 
tour  distances  produced  in  the  nearest  neighbor  phase,  and 
the  distances  of  the  final  tour  solutions  produced  by  the 
k-opt  phase  of  the  algorithm.  There  seemed  to  be  a  rela¬ 
tionship  between  the  tour  produced  by  the  nearest  neighbor 
phase  and  the  improved  tour  produce  by  the  k-opt  phase.  A 
least-squares  linear  fit  of  the  initail  tour  distance  vs. 
the  final  tour  distance  was  done  for  four  data  sets.  In 
each  case  the  fit  produced  a  positive  slope.  This  implies 
for  a  relatively  large  nearest  neighbor  tour,  the  final  tour 
from  the  k-opt  phase  will  be  relatively  large.  A  relatively 
small  nearest  neighbor  tour  will  produce  a  relatively  small 
tour  from  the  k-opt  phase.  Using  this  information  the  buoy- 
tender  algorithm  was  further  modified  to  reduce  computation 
time  without  compromising  too  much  of  its  ability  tc  produce 
optimal  solutions.  For  a  set  of  k  nodes,  k<n,  the  nearest 
neighbor  algorithm  provides  k  initial  tours.  The  set  of  k 
nodes  is  selected  at  random  from  the  population  of  n  nodes. 
Then  for  a  set  of  L  nodes,  L<k,  the  k-opt  algorithm  is  run 
to  obtain  I  final  tours.  The  set  of  L  nodes  are  drawn  from 
the  k  nodes  which  produced  the  L  shortest  initial  tours. 
The  shortest  final  tour  from  the  L  final  tours  produced  by 
the  k-opt  phase  will  be  the  algorithm^  solution  to  the 
problem.  The  final  version  of  the  algorithm  has  a  greatly 
reduced  operating  time  over  the  initial  version. 


C.  AIGOBITHH  DESCRIPTION 


In  this  section,  a  detailed  description  of  the  buoy- 
tender  algorithm  is  given.  In  this  description  references 
will  be  given  to  the  program  subroutines  so  that  the  reader 
may  associate  the  portion  of  the  algorithm  under  discussion 
with  its  operational  counterpart  in  the  program. 

Ihe  buoytender  algorithm  begins  by  reading  a  file  from 
mass  storage  containing  the  number  of  aids  (nodes)  in  the 
problem,  the  distance  matrix,  and  the  names  of  the  aids 
(subroutine  OBTAIN).  A  program  NATOTL.FOE,  which  assists 
the  user  in  creation  of  this  data,  can  be  found  after  the 
listing  of  the  buoytender  route  selection  program, 
BOUT E. FOE. 

After  obtaining  the  data  the  algorithm  selects  a  random 
sample  of  nodes,  based  on  the  size  of  n,  from  which  it  will 
create  tours  (subroutine  NODE) .  The  function  which  selects 
the  size  cf  the  random  sample  is 

Sample  Size  =  5  ♦  p2  (  Log  (  n-4  ))\  . 

The  function  is  designed  to  capture  almost  the  entire  popu¬ 
lation  for  testing  when  the  problem  is  small  and  to  capture 
only  a  small  portion  of  the  population  being  tested  when  the 
problem  is  large.  This  produces  a  sample  size  which 
increases  the  probability  of  obtaining  the  global  optimum 
when  the  problem  is  small  yet  produces  a  sample  size  which 
will  maintain  reasonable  computation  times  when  the  problem 
is  large. 

After  determining  the  sample  size,  a  random  sample  of 
nodes  is  taken.  Since  there  is  no  random  number  generator 
for  FCETEAN  installed  on  the  Coast  Suard  Standard  Terminal, 
a  pseudo-random  number  generator  was  added  to  program  (func¬ 
tion  RANDOM).  The  coding  for  the  generator  comes  from 
Wichmann  and  Hill  [Ref.  26]  ;  it  is  described  by  its 


authors  as  an  efficient  psuedo-random  number  generator 
having  a  cycle  length  of  2.78  x  10*3  .  Nodes  are  sampled  at 
random  without  replacement  from  the  integer  population  of  1 
to  n  for  use  as  starting  nodes  for  the  nearest  neighbor 
phase . 

The  algorithm  then  moves  into  the  nearest  neighbor  phase 
(subroutine  NBR)  .  A  tour  is  created  for  each  node  in  the 
random  sample,  with  the  node  from  the  sample  as  the  starting 
node  of  the  tour.  The  algorithm  selects  the  nearest  node  to 
the  starting  node  for  the  next  entry  in  the  tour  construc¬ 
tion.  The  algorithm  then  iterates  selecting  a  node  not  yet 
in  the  tour  which  is  nearest  to  the  previously  selected  node 
as  the  next  entry  in  the  tour.  This  process  continues  until 
all  nodes  have  been  selected,  at  which  time  the  last  node  is 
connected  to  the  starting  node  to  complete  the  tour,  and  the 
tour  distance  is  computed.  After  all  the  initial  tours  have 
been  constructed  they  are  ordered  by  their  total  distances 
from  shortest  to  longest.  The  five  shortest  tours  are 
retained  for  possible  improvement  in  the  k-opt  phase. 

A  decision  matrix  is  created  which  designates  which 
distance  arcs  in  the  distance  matrix  are  presently  members 
of  an  initial  tour  (subroutine  MARKD) .  The  algorithm  then 
proceeds  to  step  through  each  node  of  the  tour  testing  it 
for  k-optimality . 

The  tour  is  ’prepared*  at  each  node  to  be  checked  for 
k-optimality  (subroutine  TRPREP).  This  process  simplifies 
some  of  the  later  operations.  Here,  it  is  helpful  to 
consider  the  tour  as  a  circle  with  n  postions  on  its  circum¬ 
ference  (see  Figure  4.1),  and  to  let  positions  on  the 
circular  tour  be  denoted  by  p  (i)  where  i  =  1,-..,n.  A  node  in 
the  tour  is  represented  as  a  position  on  the  circle.  TRPREP 
places  the  starting  node  in  p ( 1 )  and  places  the  remaining 
nodes  in  the  tour  at  p(2)  through  p  (n)  where  (i)  is  the 
node’s  position  in  the  tour. 


Figure  4. 1  Circular  Representation  of  a  Tour 


Since  floating  k-cpt  algorithm  is  in  essence  a  series  of 
2-opt  iterations  a  description  of  the  siapler  2-opt  prin¬ 
ciple  Mill  fce  presented  first.  Figure  4.2  shows  the  2-opt 
principle  in  graphic  detail.  In  Figure  4.2  (A) ,  the  tour  is 
represented  in  circular  fora.  The  tour  is  now  perturbed  to 
try  and  find  an  improvement.  To  perturb  the  tour  in  the 
2-opt  procedure,  two  arcs  will  be  broken  in  the  existing 
tour  and  two  new  arcs  will  be  selected  as  incoaing  to  form  a 
new  tour.  This  provides  a  gain  function  which  can  be  used 
to  check  for  improveaent.  The  gain  function  is 

G  =  ,£(  (broken  arcs)  -  ^(incoming  arcs). 

If  G  is  positive  this  set  of  arc  exchanges  will  improve  the 
tour. 

Tc  perturb  the  tour,  one  of  the  arcs  incident  to  node 
being  checked  for  2-optimality  must  be  broken  or  removed 
from  the  tour.  In  the  buoytender  algoritha  the  node  being 
checked  for  k-optiaality  will  always  be  in  p ( 1 )  of  the 
circular  tour  and  the  arc  being  broken  will  always  be 
between  p(1)  and  p (2)  of  the  circular  tour.  An  incoming  arc 
must  then  be  selected  to  replace  the  broken  arc.  This 
incoming  arc  must  be  selected  _roa  those  arcs  with  a 
starting  node  which  is  the  same  node  as  the  end  node  of  the 


(C)  Incoming  Arc  Selected  (D)  Wrong  Incident  Arc  Broken 
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(E)  Feasible  Arc  Exchange  (F)  Circular  Tour 

Figure  4.2  Graphic  Description  of  2-opt 


previously  broken  arc.  Figure  4.2  (B)  shows  seme  of  the 
available  selections  for  the  -^coming  arc.  The  incoming  arc 


is  in  general  the  shortest  available  arc,  but  may  be  any  of 
the  available  arcs  which  maximize  G.  This  subtour  must  be 
broken  and  another  incoming  arc  must  be  selected  sc  as  to 
reconstitute  a  tour.  Therefore,  one  of  the  tour's  arcs 
incident  to  the  end  node  of  the  first  incoming  arc  must  be 
broken.  The  end  node  in  Figure  4.2  (C)  is  at  p(7).  There 
is  only  cne  of  two  possible  tour  arcs  which  may  be  broken  to 
allow  the  incoming  arc  to  become  part  of  the  tour.  Figure 
4.2  (D)  shows  the  selection  of  the  wrong  incident  tour  arc 
for  removal.  The  breaking  of  this  arc  and  the  subsequent 
reconnection  to  the  ncde  being  checked  for  2-optimality  does 
not  result  in  a  tour.  Instead,  two  subtours  are  created 
giving  an  infeasible  situation.  Figure  4.2  (E)  shows  the 
breaking  of  the  proper  incident  arc  and  subsequent  reconnec¬ 
tion  tc  the  original  node  to  create  the  'improved  tour*. 
Also  notice  that  some  of  the  arc  directions  in  the  original 
tour  must  be  reversed  to  provide  a  consistent  direction  of 
movement  through  the  tour.  Reversing  these  arc  directions 
is  the  same  as  inverting  a  sequence  of  nodes  in  the  tour. 
The  result  is  shown  in  Figure  4.2  (F)  as  a  new  'improved' 
circular  tour. 

Ihe  2-opt  procedure  is  carried  out  using  each  node  in 
the  tour  as  the  starting  node  checked  for  2-optimality.  One 
interation  of  the  2-opt  procedure  through  all  n  nodes  will 
produce  a  tour  which  is  2-optimal. 

The  k-opt  algorithm  is  very  similar  to  the  2-opt  proce¬ 
dure  described  above.  Instead  of  just  reconnecting  the  tour 
as  the  2-opt  procedure  does  (p(6)  to  p  ( 1  > )  ,  the  k-opt  checks 
to  see  if  another  incoming  arc,  starting  at  the  end  ncde  of 
the  last  broken  arc,  can  be  selected  to  further  improve  the 
tour  (Figure  4.3).  The  procedure  interates  until  no  further 
improvement  can  be  found.  A  k-opt,  using  an  example  with  k 
equal  to  four,  might  look  like  Figure  4.4  The  floating  k-opt 
algorithm  has  the  ability  to  vary  k  during  the  k-opt 


Figure  4.3  K-opt  Selecting  Next  Incoming  ire 


Figure  4.4  Example  of  4-opt 


solution  process.  This  is  sore  desirable  than  fixing  k  in 
the  solution  process.  For  example,  if  k  vere  fixed  at  four, 
the  algorithm  can  only  check  four  arc  exchanges  for  improve- 
sent  but  the  problem  may  need  a  five  arc  exchange  to  produce 
the  optimal  solution.  By  allowing  the  k-opt  to  vary  k,  the 
probability  of  the  algorithm  finding  the  optimal  solution  is 
increased.  The  floating  k-opt  is  bounded  between 
2-optimality  and  n-optimali ty.  The  k-opt  cannot  be  greater 
than  n-optimal  because  an  n-opt  exchange  creates  the  optimal 


The  k-opt  procedure  in  the  buoytender  algorithm  operates 
by  checking  each  node  in  a  tour  for  k-optimality.  The  algo¬ 
rithm,  using  subroutine  TR PREP,  creates  the  circular  tour 
for  the  node  being  checked  for  k-optimality.  The  circular 
tour  is  prepared  so  that  the  initial  arc  being  broken  will 
always  be  between  p(1)  and  p(2).  The  algorithm  checks  the 
circular  tour  at  the  starting  node  for  k-optimality.  If  no 
improvement  can  be  made  breaking  the  first  arc  between  p ( 1 ) 
and  p  (2)  then  the  circular  tour's  direction  is  reversed 
(subroutine  REVTOR)  so  that  the  alternate  incident  arc  is 
placed  between  p  (1 )  and  p(2)  of  the  circular  tour  and  this 
tour  is  checked  for  k-optimality. 

The  breaking  of  the  arc  between  p(1)  and  p  ( 2)  requires 
an  incoming  arc  to  replace  it  (subroutine  SELCTT)  .  The 
algorithm  selects  the  five  shortest  available  arcs  which 
originate  at  position  two  of  the  circular  tour  and  terminate 
at  another  position  on  the  circular  tour.  The  arc  from  p  (2) 
to  p(3)  is  not  available  since  it  is  already  in  the  tour, 
and  the  arc  from  p(2)  to  p  (1)  is  not  available  since  it  is 
the  arc  just  broken  to  improve  the  tour.  For  each  of  the 
selected  incoming  arcs  another  arc  must  be  broken  to  elimi¬ 
nate  a  subtour  (as  in  the  2-opt  procedure)  .  It  is  then 
necessary  to  reconnect  the  arcs  to  reconstitute  the  tour. 
This  results,  for  each  of  the  selected  arcs,  a  set  of  two 
incoming  arcs  and  two  outgoing  arcs.  Using  the  previously 
defined  G,  each  of  the  five  sets  of  arcs  is  checked  to  see 
which  maximizes  G.  If  G>0  then  the  exchange  or  inversion  is 
made  (subroutine  ADJTCR)  ,  otherwise  no  exchange  is  made. 

If  an  exchange  is  made  the  algorithm  goes  back  and 
searches  for  a  new  incoming  arc  from  the  end  point  of  the 
last  broken  arc.  Subroutine  ADJTOR  has  inverted  the  tour  so 
that  the  last  arc  used  to  reconnect  the  tour  is  now  between 
p ( 1 )  and  p  (2)  of  the  circular  tour.  Subroutine  SELCTY  is 
again  used  to  select  the  next  best  available  incoming  arc  as 


described  above.  The  algorithm  interates  through  selecting 
an  incoming  arc  and  inverting  the  tour  as  long  as  improve¬ 
ment,  i.e.  positive  gain,  can  be  found. 

When  improvement  can  no  longer  be  found  the  direction  of 
the  tour  is  reversed  at  the  last  node  where  improvement  was 
found,  and  the  tour  is  tested  again  for  possible  further 
improvement.  When  nc  further  improvement  can  be  found  and 
the  reversed  tour  has  also  been  tested,  the  algorithm  then 
increments  to  the  next  node  in  the  tour  and  starts  the  k-opt 
procedure  again  testing  this  node  for  k-opt imality.  Each 
node  in  the  tour  is  used  as  a  starting  node  in  the  k-opt 
procedure.  When  all  nodes  in  the  tour  have  been  tested  for 
k-optimality  the  k-opt  phase  ends  and  the  resulting  tour  is 
stored. 


After  each  of  the  five  shortest  initial  tours  have  been 
tested  for  k-optimality,  the  shortest  of  the  improved  k-opt 
tours  is  selected  as  the  solution  to  the  buoy tender  problem. 
Figure  4.5  is  a  flowchart  of  the  buoytender  algorithm. 

The  next  chapter  discusses  the  results  of  some  test 
problems  solved  by  the  algorithm  and  makes  concluding 
remarks  about  the  algorithm. 


Figure  4.6  Algorithm  Flowchart  Section  2 


V.  PROGRAM  EESOLTS  AND  CONCLUSIONS 


This  chapter  will  discuss  some  results  obtained  from 
using  the  tuoytender  route  selection  program,  including 
results  from  an  actual  problem  from  the  U.S.  Coast  Guard 
Cutter  BLACKHAW.  All  program  results  were  obtained  using 
the  Coast  Guard  C3  Standard  Terminal.  The  chapter  will 
conclude  with  some  observations  regarding  the  buoy  tender 
problem,  the  algorithm,  and  the  test  results. 

A.  PBOGEAH  EESOLTS 

To  program  and  test  the  buoy tender  algorithm  it  was 
necessary  to  obtain  problems  with  known  optimal  solutions. 
Initially,  a  very  simple  and  small  problem  with  nine  nodes 
was  created.  This  test  problem’s  optimum  tour  is  the  most 
obvious  tour  obtained  by  inspection  of  the  problem  since  the 
nodes  are  arranged  in  an  almost  circular  pattern. 
Connecting  these  nodes  following  this  nearly  circular  path 
results  in  the  optimum  solution.  Several  alternate  tours 
were  checked  by  hand  computation  and  were  found  to  be  subop- 
timal.  Therefore,  the  near  circular  path  for  this  test 
problem  is  assumed  to  be  optimal.  When  the  program  was 
tested  with  this  test  problem,  the  program  produced  this 
hypothesised  optimum  solution  with  a  computation  time  of  13 
seconds. 

With  the  knowledge  that  the  program  appeared  to  be  oper¬ 
ating  correctly  on  tie  initial  test  problem,  further  testing 
was  necessary  to  check  the  program's  ability  to  solve  other 
problems.  Since  the  tuoytender  route  selection  problem  is  a 
symmetric  case  of  the  TSP,  the  literature  was  searched  for 
test  problems.  Three  articles  [Ref.  5,20,27]  yielded  eight 


symmetric  TSP's  of  various  sizes.  These  eight  literature 
problems  ranged  in  size  from  5  to  57  nodes  with  known 
optimal  answers.  The  program,  when  tested  with  these  prob¬ 
lems,  produced  solutions  which  ranged  from  the  optimal  to 
within  3.5  percent  of  the  optimal  distance. 

Another  problem  was  created  by  the  random  selection  of 
70  pcints  on  a  Euclidean  plane.  A  distance  matrix  was 
computed  for  these  70  nodes  and  this  problem  was  solved  by 
the  program.  The  problem's  optimal  solution  is  not  known, 
but  the  problem  was  run  to  ensure  the  program  could  handle 
the  maximum  problem  for  which  it  was  designed,  and  to 
produce  a  sample  solution  time  for  this  maximum  sized 
problem. 

With  evidence  that  the  program  was  solving  test  prob¬ 
lems,  a  further  test  was  conducted  using  an  actual  buoy- 
tender  routing  problem,  with  data  provided  by  the  U.S.  Coast 
Guard  Cutter  BLACKHAW  (WLB-390)  ,  homeported  in  San 
Fransisco.  The  aids  to  navigation  run  data  was  for  a  trip 
scheduled  from  9  April  1984  to  16  April  1984.  With  the 
experience  and  historical  data  this  vessel  has  available 
from  servicing  these  aids  for  many  years,  it  is  hypothesized 
that  the  route  which  was  scheduled  for  this  ATON  run  would 
be  optimal  or  nearly  optimal.  The  planned  route  for  the 
ATON  run  had  a  distance  of  455  nautical  miles.  When  this 
problem  was  entered  and  solved  by  the  buoytender  program, 
the  solution  was  also  455  nautical  miles.  Although  the 
distances  were  the  same,  the  scheduled  route  and  the  program 
route  differed  slightly  in  their  sequence  for  visiting  the 
aids.  This  difference  may  be  attributed  to  the  fact  that, 
as  stated  in  Chapter  II,  the  triangle  inequality  may  be 
required  to  be  a  strict  equality  on  some  of  the  arcs.  This 
difference  is  possible  since  the  buoytender  may,  due  to  the 
strict  equality  on  some  arcs,  be  required  to  'pass'  a  previ¬ 
ously  serviced  aid  to  complete  the  route.  Figure  5.1  shows 
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a  five 

node  example  of 

this  situation 

.  This  suggests  there 

say  be 

several  sequences  of  nodes 

which  will  provide  an 

optimal 

solution  to  the 

buoytender  problem. 

TABLE  1 

Tours  : 

for  BLiCKHAH  Problem 

Route  as  Scheduled 

Route  Produced 

hi 

the  Buoy  tender 

by  the  Proqram 

1 

YB  ISLAND 

1 

YB  ISLAND 

2 

MONTAR A  LNB10A 

20 

MOSS  LBB  MLA 

3 

AN  ISLAND  LWE8 

19 

MONTRY  BY  LB  B 

4 

SANTA  CROZ  ME 

18 

MONTRY  HBR  MB 

5 

P  ELANCAS  LWE4A 

17 

PT  PINOS  LHB2 

6 

SAN  SIMEON  LEB1 

16 

PT  CYPRESS  LGB6 

7 

VON  HELM  R  LGE4 

5 

P  BLANCAS  L WB4A 

8 

ESTERO  BY  GB10E 

6 

SAN  SIMEON  LBB  1 

9 

SOOZA  R  LG  B 

7 

VON  HELM  R  LGB4 

10 

NESTDAHL  R  LEE1 

15 

ESTERO  BY  LHB 

11 

£T  SAN  LOIS  LNE3 

8 

ESTERO  GB10E 

12 

LANSING  R  B 

14 

MORRO  BY  LBB1 

13 

FT  BOCHON  LSE2 

13 

PT  BOCHON  LHB2 

14 

MORRO  BY  LBB1 

10 

NESTDAHL  R  LBB1 

15 

ESTERO  BY  LHE 

9 

SOOZA  R  LGB 

16 

FT  CYPRESS  LGE6 

12 

LANSING  R  B 

17 

FT  FINOS  LHB2 

11 

PT  SANLOIS  R  B 

18 

MONIRY  HBR  ME 

4 

SANTA  CROZ  MB 

19 

MONTRY  BY  LB  E 

3 

AN  ISLAND  LWB8 

20 

MOSS  LBB  MLA 

2 

MONTARA  LNB10A 

1 

YB  ISLAND 

1 

YB  ISLAND 

Figure  5. 1  Multiple  Optiial  Solutions 


TABLE  2 

Results  of  Test  Problems 


Problem  6 
(aef.) 

Size 

Solution 

Time  (Sec) 

Tour  Dist 
££oq ram 

:.  (Miles) 
Optimal 

Solution 

Quality* 

Xarg(27) 

Tes t  (auth. ) 

5 

5 

148 

148 

100* 

9 

13 

22.95 

22.95 

100* 

Barachet  (27; 

1  10 

17 

378 

378 

100* 

3LACKHA  W 

20 

50 

455 

455 

** 

Croes  (20) 

20 

52 

254 

246 

103.  25* 

HeldKarp  (5) 
Kar  g  (27) 

25 

72 

1711 

1711 

100* 

33 

121 

10929 

10861 

100.  62* 

Dant2ig  (27) 
HeldKarp  (5) 

42 

188 

705 

699 

100.  86* 

48 

238 

11847 

1 1470 

103.28* 

Kara  (27) 

57 

340 

13367 

12956 

103.  17* 

Max  (auth. ) 

70 

494 

6791. 07 

unk. 

.  unk 

*  100*  is  optimal; 

103.28*  is  3 

.28*  greater  than 

optimal. 

**  The  distance  is 

hypothesized 

to  be  optimal. 

Figure  5.2  Least  Squares  Fit  of  Size  vs.  Solution  Tine 

Table  2  presents  a  summary  ox  the  eleven  test  problems 
shoving  problem  sizes,  solution  times,  and  the  quality  of 
the  solutions  produced  by  the  program.  All  of  the  run  times 
are  reasonable  by  the  criteria  stated  in  Chapter  II.  The 
largest  problem,  70  nodes,  had  a  solution  time  of  494 
seconds,  or  roughly  a  little  more  than  eight  minutes. 


Computation  times  for  the  eleven  problems  were  plotted 
against  problem  size,  as  shown  in  Figure  5.2  A  least-squares 
polynomial  fit  of  the  data  points  is 

Time  =  .584  ♦  .62n  ♦  .0914n2  ♦  .00000994n* 
and  this  curve  is  also  shown  in  the  figure. 


E.  CCNC1USIOHS 

The  buoytender  route  selection  algorithm  provides 
another  useful  way  of  scheduling  ATON  routes.  As  a  tool  for 
the  operator  to  assist  in  planning  the  route  to  be  taken, 
the  program  should  produce  optimal  or  near  optimal  solutions 
for  prcblems  up  to  a  size  of  70  aids.  The  user  may  accept 
this  route  as  is  or  may  desire  to  modify  the  route  based  on 
circumstances  requiring  human  judgement.  Appendix  A  is  a 
user's  guide  for  the  matrix  utility  and  buoytender  route 
selection  programs. 

The  buoytender  route  selection  algorithm,  with  the 
nearest  neighbor  and  k-opt  heuristics,  provides  a  quick  and 
satisfactory  solution  to  a  symmetric  traveling  salesman 
problem  cn  problems  up  to  a  size  of  70  nodes.  The  problem 
may  have  several  routes  with  different  sequences  for 
visiting  the  aids,  tut  the  routes  may  all  have  the  same 
distance. 

It  is  hoped  that  this  algorithm  will  be  of  use  to  the 
fleet  cf  O.S.  Coast  Guard  buoy  tenders.  As  a  decision  aid  it 
may  help  the  operators  of  buoytenders  obtain  optimal  or 
nearly  optimal  routes  to  service  their  aids,  possibly  saving 
time  and  fuel. 


IPPEHDIX  A 

USEE’S  60 IDE  TO  PROGRAMS 


The  intention  of  this  appendix  is  to  provide  instruc¬ 
tions  for  the  use  of  the  buoytender  route  selection  program 
and  utility  programs.  It  should  be  noted  that  prior  to 
using  the  buoytender  route  selection  program,  the  user  will 
need  to  use  the  matrix  utilities  program  listed  at  the  end 
of  this  appendix  to  prepare  a  data  file. 

The  buoytender  route  selection  program  and  the  matrix 
utilities  program  are  written  in  Convergent  FORTRAN. 
Convergent  FORTRAN  is  FORTRAN  77  compatible.  The  programs 
should  be  typed  in  as  presented  and  linked  to  the  operating 
system.  The  Convergent  FORTRAN  manual  has  instructions 
concerning  compiling  and  linking.  Recommended  names  for  the 
run  files  are  ROUTE. BCN  and  MATUTL.RUN.  Both  programs  will 
operate  cn  either  a  INS  or  AWS  station  with  264 K  memory  or 
more.  The  programs  run  about  three  time  slower  on  an  AWS 
terminal  than  on  an  I  '  terminal. 

The  user  will  be  required  to  enter  text  and  numerics 
into  the  programs.  All  numeric  entries  are  to  be  integer 
entries  (no  decimal  point)  with  the  exception  of  the 
distance  entries.  All  distance  entries  require  a  decimal 
point  entry  for  proper  input. 

MATUTL.RUN  is  capable  of  entering  problems  up  to  a  size 
of  100  aids.  ROUTE. RON  is  designed  to  handle  problems  up  to 
a  size  of  70  aids.  MATUTL.RUN  allows  the  entry  of  a  large 
matrix  of  aids  and  includes  a  utility  which  allows  the  user 
to  select  a  subset  of  this  large  matrix  for  use  in  the 
program  ROUTE. RUN.  This  allows  the  user  to  type  in  one 
large  matrix  and  then  create  smaller  matrices  from  the 
larger  as  needed.  This  feature  wi.l  preclude  the  user  from 


typing  in  a  distance  matrix  every  time  he  desires  to  run  a 
problem. 


A.  USE  CF  NATUT1.  BOH 

The  matrix  utilities  program  is  designed  for  creating 
and  editing  data  files  which  will  be  used  in  the  buoytender 
route  selection  program.  A  data  file  will  contain  the  size 
of  the  problem  n,  the  distance  matrix  for  the  problem,  and 
the  names  of  the  aids  to  be  visited. 

The  size  of  the  problem,  n,  will  be  the  number  of  aids 
to  be  visited  plus  one.  The  plus  one  is  for  the  port  from 
which  the  buoytender  will  start  and  end  the  ATON  run. 

To  create  the  distance  matrix  some  extensive  plotting  is 
required.  Most  buojtenders  have  a  set  of  tracklines  they 
regularly  follow  to  visit  their  aids.  In  most  cases  these 
existing  tracklines  will  be  sufficient  to  use  in  the 
distance  matrix.  To  create  a  distance  matrix  for  entry  into 
a  data  file  it  is  recommended  the  user  use  the  format  shown 
in  Figure  A-1  .  Distances  should  be  computed  in  nautical 
miles  and  may  range  from  .0  1  nm  to  9999.99  did.  Figure  A.  1 
is  an  example  of  a  problem  of  size  nine. 

Since  the  problem  is  symmetric,  the  user  need  only 
calculate  and  enter  the  distances  as  shown  in  Figure  A.  1  . 
The  matrix  utilities  will  automatically  complete  the  rest  of 
the  matrix  entries.  The  distance  you  should  compute  should 
be  the  shortest  possible  distance  from  the  row  entry  to  the 
column  entry.  Using  the  example  in  Figure  A.1  the  entry  in 
row  3,  column  6  is  the  shortest  possible  trackline  from  aid 
A  to  aid  E.  Once  the  entries  for  the  distance  matrix  have 
been  prepared  the  user  is  ready  to  use  MAT UTL .RUN . 

M ATUTL. RUN  has  seme  simple  menus  and  some  descriptive 
comments  to  help  the  user.  The  menus  require  a  numeric 
response  for  the  selection  of  an  option.  The  first  menu 
presents  the  utilities  available  to  the  user.  Figure  A. 2  . 
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Figure  A.1  Sample  Distance  Matrix 
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1  -  create  a  distance  matrix 

2  -  extract  a  small  matrix  from  a  larger  matrix 

3  -  correct  a  previously  stored  matrix 

4  -  create  a  hardccpy  of  a  distance  matrix 

5  -  list  names  with  their  row  number  in  the  matrix 

6  -  exit  this  pregram 

What  is  your  selection? 


Figure  A. 2  MATUTL.EUN  Main  Menu 

1.  Create  a  Distance  Matrix 

This  utility  is  for  the  initial  entry  of  a  distance 
matrix  into  a  data  file.  The  user  will  be  asked  for  the 
size  of  the  problem,  the  homeport  and  aids  names,  and  the 
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distances.  The  size  of  the  problem  is  the  number  of  aids 
plus  one,  for  homeport.  The  utility  will  next  ask  for  names 
of  homeport  and  the  aids  to  be  visited.  These  names  will  be 
associated  with  a  row  and  column  number  of  the  distance 
matrix.  The  names  are  text  entries  up  to  a  maximum  of  15 
characters. 

The  utility  will  next  ask  for  distances.  It  will 
ask  for  the  distance  from  a  row  number  to  a  column  number. 
If  the  user  sets  up  the  distance  matrix  as  shown  in  Figure 
A. 1  he  will  be  able  to  then  just  read  each  row  across  its 
columns  while  entering  the  data.  These  distances  must  be 
real  numbers.  The  distances  should  be  nautical  miles  and 
may  range  from  .01  nm  to  9999.99  nm.  Since  these  entries 
are  real  numbers  there  must  be  a  decimal  point  in  each 
entry.  If,  for  example,  the  user  needs  to  enter  two 
nautical  miles,  he  should  enter  it  as  2.  <return>.  Be 
careful  with  these  entries;  if  a  decimal  point  is  not 
entered  the  distance  will  be  off  by  a  factor  of  100.  For 
quickly  and  easily  entering  data  the  numeric  keypad  on  the 
keyboard  may  be  used  and  <next>  may  be  used  in  place  of 
<return>  to  enter  the  data.  If  an  error  is  made  while 
typing  in  numbers  and  <return>  or  <next>  have  not  yet  been 
pressed,  then  the  error  may  be  corrected  by  using  the  back¬ 
space  key  and  retyping  the  number.  If  an  error  is  made  and 
<return>  or  <next>  have  already  been  pressed,  then  the  user 
should  note  the  row  and  column  number  of  the  error  and  use 
the  correction  utility  to  correct  the  entry. 

When  the  user  finishes  typing  all  entries  the 
program  will  prompt  for  a  filename  in  which  the  data  will  be 
stored.  This  name  may  be  up  to  15  characters  long.  If  for 
some  reason  the  user  needs  to  end  the  terminal  session 
before  the  entire  matrix  is  entered,  data  entry  may  be 
stopped  by  typing  in  a  negative  distance.  The  user  will 
then  be  asked  for  a  filename  for  the  partial  data  file.  To 
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complete  the  partial  data  file  the  user  must  use  the  correc¬ 
tion  utility. 


2.  Extract  a  Matrix 

This  utility  is  designed  so  tnat  the  user  may  choose 
a  subset  of  aids  from  a  larger  matrix.  The  utility  will 
create  the  proper  data  file  for  a  subset  of  aids. 

The  user  should  have  an  idea  of  all  possible  aids 
which  may  be  visited  on  a  general  buoytender  run  (i.e.  the 
spring  run,  the  north  run,  the  inside  run,  etc.) .  The  user 
should  create  a  distance  matrix  using  the  create  utility  of 
all  aids  which  could  possibly  be  visited  on  a  general  run. 
Usually  cnly  some  portion  of  these  aids  will  be  visited  on 
an  actual  run  (for  example  on  a  particular  spring  run  only 
90*  of  the  aids  which  might  be  visited  on  this  run  need  to 
be  visited) .  In  this  case  the  user  creates  a  matrix  of  all 
aids  which  might  be  visited  on  a  run  and  then  each  time  this 
run  is  scheduled  the  user  can  then  use  the  extract  utility 
to  quickly  and  easily  create  a  data  file  for  those  aids 
which  need  to  be  visited. 

Before  using  the  ex  tract  utility  it  is  suggested  the 
user  first  use  the  list  utility  on  the  larger  or  initial 
data  file.  This  will  provide  the  user  a  list  of  aid  names 
in  the  data  file  and  present  them  with  their  corresponding 
row  number. 

The  extract  utility  will  prompt  the  user  for  the 
name  of  the  data  file  containing  the  general  run*s  aids. 
The  user  is  then  asked  for  the  size  of  the  new  problem 
whereby  the  user  enters  some  m,  m  less  than  n  in  the  initial 
data  file.  The  utility  then  asks  the  user  m  times  for  row 
numbers  from  the  initial  matrix.  After  m  row  numbers  have 
been  entered  the  utility  creates  the  new  data  file  and  askes 
the  user  for  a  filename  for  this  new  data  file. 


3.  Correct  a  Previously  Stored  Matrix 

This  utility  is  to  correct  entries  in  a  distance 
matrix  or  tc  complete  a  partially  entered  data  file.  The 
utility  has  a  menu  which  askes  if  the  user  wishes  tc  make 
single  entry  or  seguential  entry  corrections.  The  single 
entry  mode  asks  the  user  for  the  row  number  and  column 
number  of  the  entry  desired  to  be  corrected.  The  user  then 
is  asked  for  the  new  distance  and  is  returned  to  the  correc¬ 
tion  menu  where  he  may  select  to  make  another  correction  or 
exit  the  utility. 

The  seguential  corrections  are  similar  to  the  entry 
of  distances  in  the  create  utility.  The  user  is  asked  for 
the  starting  row  number  and  starting  column  number  where 
corrections  are  to  begin.  The  utility  asks  for  the 
remaining  column  entries  on  the  starting  row  and  then  will 
proceed  to  the  next  row  until  entries  are  terminated.  The 
utility  automatically  terminates  at  the  end  of  the  matrix  or 
the  user  may  terminate  the  seguential  mode  early  by  entering 
a  negative  number.  If  a  negative  number  is  entered  it  will 
not  appear  in  the  distance  matix.  Terminating  the  seguen¬ 
tial  mode  returns  the  user  to  the  corrections  menu  for 
another  correction  or  to  exit  the  correction  utility. 

Upon  exiting  the  correction  utility  the  user  is 
asked  for  a  filename  for  the  corrected  file.  This  filename 
may  be  up  to  15  characters  long. 

4 .  Fre par e  a  Distance  Matrix  for  Hardcopy 

This  utility  is  for  preparing  a  hardcopy  of  a 
distance  matrix  from  a  stored  data  file.  The  data  files  are 
saved  as  seguential  files,  (i.e. ,  the  data  files  are  a 
string  of  numbers) .  This  utility  will  construct  a  file 
which  will  resemble  a  matrix.  Due  to  the  width  and  length 
limitations  of  the  printer  and  paper,  the  matrix  is  printed 


in  sections  which  will  have  to  be  cut  and  pasted  together  to 
display  the  whole  matrix.  A  hard  copy  of  the  whole  matrix 
will  provide  the  user  with  a  means  of  checking  to  see  if 
data  entries  are  correct.  If  an  entry  error  is  found,  the 
matrix  also  helps  the  user  to  locate  the  row  and  column 
number  of  the  error. 

The  utility  asks  the  user  for  the  name  of  the  data 
file  from  which  a  distance  matrix  is  needed.  The  utility 
creates  the  hardcopy  file  and  then  prompts  the  user  to 
obtain  a  hardcopy  by  using  the  FORMAT  command  to  print  out 
file  MATEET .DATA . 

5.  list  Aid  Names 

This  utility  provides  a  hardcopy  of  the  aid  names 
from  a  data  file.  This  may  be  useful  to  the  user  when  using 
the  extract  utility.  The  user  is  asked  for  the  name  of  the 
data  file  from  which  the  list  shall  be  made.  The  utility 
then  prepares  a  file  containing  the  row  numbers  and  their 
associated  aid  names.  The  user  is  prompted  to  use  the 
FORMAT  command  to  obtain  a  printout  of  LIST. DATA. 

These  matrix  utilities  are  for  assisting  the  user  in 
creating  and  manipulating  data  files  and  distance  matrices. 
A  correct  data  file  is  necessary  for  the  proper  operation  of 
the  buoytender  route  selection  program,  ROUTE. RUN. 

B.  USE  CF  ROUTE. RUH 

The  user  of  ROUTE. EUN  is  very  simple  once  a  data  file 
has  been  created.  The  user  runs  the  ROUTE. RUN  program  and 
he  is  prompted  for  the  name  of  a  data  file.  The  user  is 
then  asked  for  the  name  of  an  output  file  for  the  results. 
The  program  obtains  the  size  of  the  problem  from  the  data 
file  and  will  compute  an  approximate  solution  time.  This 
approximate  solution  time  is  based  on  results  from  some  test 
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problems  solved  by  the  program  and  it  assumes  the  user  has 
version  8.01  of  the  operating  system.  If  another  version  of 
the  operating  system  is  used  the  solution  time  may  differ. 

After  displaying  the  approximate  solution  time  the 
program  reads  the  rest  of  the  data  file  and  proceeds  to  do 
the  necessary  computations.  The  user,  depending  on  the  size 
of  the  problem,  will  have  just  enough  time  to  get  a  cup  of 
coffee  before  obtaining  the  solution.  The  solution  will  be 
displayed  on  the  screen  and  will  also  be  placed  in  the  file 
previously  named  for  results,  from  which  a  hardcopy  may  be 
made.  The  user  now  has  a  possible  route  to  take  for  that 
particular  aids  to  navigation  run. 


&PPENDII  B 

SABP1E  OOIPOT  OP  EOOYTENDEB  BOOTE  SELECTION  PBOGBAM 

The  solution  for  data  file  3 LACKHA W .  M  AT 
Total  tour  distance  =  455.00 
The  tour  is  as  follows 

1  YB  ISLAND 

20  MOSS  LBB  MLA 
19  MONTRY  BY  LB  3 
18  MONTEY  HBR  MB 
17  PT  PINOS  LBB2 
16  PT  CYPRESS  LGB6 

5  P  BLANCAS  LWB4A 

6  SAN  SIMEON  L3B1 

7  VON  HELM  R  LG34 

15  ESTERO  BY  LW  B 

8  ESTERO  GB10E 
14  MORRO  EY  LBB  1 
13  PT  3UCHON  LW 32 

10  KESTDAHL  R  L3B1 

9  SOOZA  R  LG B 
12  LANSING  R  B 

11  PT  SANLUIS  R  3 

4  SANTA  CRUZ  M3 

3  AN  ISLAND  LWB8 

2  MONTARA  LWB10A 
1  YB  ISLAND 


8 


129 


132 


112 


1  an 


143 
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APPENDIX  C 

LISTING  OF  BOOXTENDEB  BOOTE  SELECTION  PROGRAM 

MAIN  PROGF^'"1  ***  PUOYTENOF.P  ROUTE  SFLFCTI7'!  *** 
IMPLICIT  CHAPACTFP* 1 5  (A) 

REAL  G,TDIF,TCRDIS,DIF (15) 

INTEGER  I,J,L,C,P,  MUM ,  NO  D ,  F  AN  ( 1  5) 

C  HA  R ACT  E  P*  1  5  F  N7-  '■*  E  ,  CNA  v  E 
C CM.MCN  /NAME/  A NAT F ( 7 2 ) 

C  ON MON  /W  A L Y /  N ,  I  TOUR (71) 

COMMON  /  TDST/  TO  (7  3, 7  3) 

COMMON  ,/INVT/  IT  (14^) 

COMMON  /CIST/  0(73,72) 

COMMON  /FAMT/  IF TOR (15, 71) 

COMMON  /Tv TP/  ITPTCP  (14?) 

CALL  03 T A  I N ( F NAM E , ON AM E , M ) 

0  PEN  (  3  ,  F  I  LE=C  V,AM  E ,  STATUS  =’ N  E'.'.”  ,  FCRM=  '  F ORMATTEO 1  ) 
TORDIS=9Q99Q . 99 
CALL  NODE (N , NUM , FAM) 

CALL  N’C R  ( N  UM  ,  SAM  ,  E  I  F  ) 

DO  7  L  =  1 , 5 

DC  8  J  =  1 ,  H  + 1 

I  TOUR ( J ) = I STOP. ( L , J ) 

CONTINUE 
CALL  MAP ED ( ) 
nfl  1  a  a  n  _  1  M 

CALL  TRPF  EP  f C) 

R  =  3 

1=1 

IT ( 1 ) = I  TOUR ( 1 ) 

IT (2) = I TOUP (2) 

:-LL  FELCTY ( I , C,G, NOD) 

1=1  +  1 

IF(G.Lr.  (  3 . 2  3  9 1  )  )  G  OT  0  1 1  01 
CALL  ADJTCR ( I , P, NOD) 

GOT°  133 

IF ( ? . EC. 1)GCT0103 
CALL  DFVTCP ( R ) 

GOTO  1 2 
T  D I  S  =  9  . 9 
DO  14  2  J  =  1 , M 

TDI S=TD I F  +D (I  TOUR (J)  , I  TOUR (J  +  l)  ) 

C  ONT I v OF 

IF (TOIF . GF. T^poi ? ) GOTO? 

T''?DIF=7DIS 
DO  5  J  =  1  ,  N  j-  1 

I  STOP  ()  ,  J)  =1  TOUR  (J) 

Ts  S  *  ^  +  »!MC“ 

Cr-MT  I  r?JE 
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~C  23  1=1,  »l 

ITCCR  ( I  )  «  I  s  "OR  ( 1  ,  I) 

ccmt  :  nve 

CALL  TPPREPf!) 

WPITE  (*  ,2?)  F':A«F 
WRITF  (3 , 2  A)  FN’ANE 
W  R I T  F  (  *  ,  3  •? )  T  0  F  .0 1 F 
W  P  1 7  E  ( 3  ,  3  ? )  TOPCIS 
WRITE (*,?!) 

W FITE (3,31) 

WRITE (*, 32)  (  TTCCR  ( I  )  ,  A  *’A'*F  ( I  TOL’F  (  T  )  )  ,  I  =  1  ,  v’+ 1  ) 
WRITE  (3  ,  32)  (ITC-yo  ( I  )  ,  VIA'1  E  ( I  T^MR  (  T  )  )  , 1=1 , N+l ) 
CLOSE  (3  ) 

STOP 

-  OPXAT  ( 1  X ,  ’  The  solution  for  Data  file  1  ,  AO 
FCRy  AT  ( 1 X ,  ’  To  t  1 1  tour  'distance  =  ',FE.2) 

FORM-'  T  ( 1  v,  '  Th  a  tour  is  as  follows:') 

EOP.VAT  (IX,  I  3 ,  5 X  ,  A ) 

END 

S  'JP  P  OUT  I V  E  CP  T  A  I N  (  Fv '  v  F  ,  0\'AM  E , N ) 

IN  FLIC  IT  CHARACTER* 1 5  (A) 

INTEGER  >1 
REAL  TIMF 

CH& F  A C T E P *  1  "  " \A v E ,  0 v c 
CC.y.ycM  ,/;;m*e/  ava**e  i~n ) 
ccyy cn  /t t st/  r(?n,7^) 

W  P  IT  E  (  *  ,  6  ) 

P  E A  E  (  *  ,  7  )  E'7-.v  F 

WRITE  {■*•,=») 

PEAD(*,7)  rvAvp 

OPEN"  (4  ,  E  I  L  r  =  E  *' a  ”  E  ,  ~  o  p  v  =  •pf^v^TC'  •  ^ 

PEAT  (4  ,  'A,  E p p  =  1 ,1 )  *; 

TIME=.3n4  +  .o2*N-*-.Rc)l^*'!**2+.ni‘'1Or,',Rn4*N**~ 
'•RITE  (*,!'’)  T  Iv  E 
’WRITE  (  *  ,  i  1) 

PEAT  (4  ,  12,  EPF*  ■  A)  (  ( 0  ( T  ,  J )  ,  J  =  1 ,  “ ' )  ,1  =  )  ,N) 

F  PAD (4 , 7, TFP  =  L  M  (A':AVF(  I  1  ,  1  =  1 ,  N) 

CLOSE  (4  ) 

77  T~  T  •  r-> 

■  RITE  (*  ,  1-) 

C  LC  S  V  ( 4  > 


E  ’’"“AT  (IF,'  'An  ~  or  ■..■Wat  f  i  1  -•>  ir;  is  f.v»  3  i  s 
FT PYAT  (  A) 


POE  *  A  T ( 1 y  '"l! 


for  the  r  ~sul ts?  ' x ) 


m  ro 


VKV  f!'  =  TD  (  I  STOP  ( L  ,  I )  ,  I  TEMP) 

DC  2 2  J=2,N 

lb'  (TC  ( I  STOP  (L  ,  I)  ,  J)  .GO. TEMPO) GOTO 2 2 
IT  DMP=J 

TE''FO=TD  (I  STOP  (L,  I)  ,  J) 

CONTINUE 

I  STOP  ( L  ,  I  +1 )  =  I  TF*'  P 
D I S ( L ) =0  I S ( L ) +  TFMPO 
CO  23  K = 1 , N 

TO  (K  ,  ITF'-P)  aa 

CONTINUE 
C  O  N  T  I  N  U  F 

ISTCRl(Lf  N+l)=ISTOR  (L,  l) 

CIS (L) =DIS (L) +P ( I  STOP  (L, M)  , I  STOP  (L,  I) ) 
CONTINUE 
N  0  F  X  =  1 

DO  210  1=1,  MUM - 1 

I  F  (  C  I  S  ( I  )  .FP.DIS(I-t-l)  )DIS(l+l)=?n9°?.r,9 
IF(DIS(I) .Lr.CTS(I+l) ) GCT02 1 a 
N  C  E  X  =  0 

T  CM  PD  =  D  !  S  ( I A 1 ) 

CIS (I+1)=DIS  (I) 

CIS ( I ) = TEMPO 
DC  22  ^  J=1,N+1 

I  TF‘'P=I  STOP  (  T  +1 ,  J) 

I  STC  P ( I  +  1 , J ) = I  STOP (I ,J) 

ISTOp  ( I  ,  J)  =  I  T  F -v1  P 

*  ""p 

Z  I  S’T’P 

IF(NOEX.  NF.  1 )  GOTO  2  5 
F  ETL’F  M 


SUBROUTINE  TRPREP(C) 

INTEGER  C ,  J ,  K 
COMMON  /MOL V./  N  ,  ITCUF  (71) 
COMMON  /TV  TP/  ITPTPPfl^'1) 
J=J 

DC  4'1  F=1,N 

ITPTCF (K ) =T  TCf«  (K) 
ITPT^P (K+N)  =  I  TOUR  ( Y ) 

IF  (  I  TPTOR  (K )  .  F'2.  C)  .7  =K 

''••rf-  TrfT 

dc  n  v 

I  TOUR  (?')  =  ITPTr0  (  J*”-l  ) 


I  M  'i  i  r^/  *J*  ■!*  ■.^  V  wm-'  ^.i  -sr  -v  -w  -  *j**zt*  ’■-"' w.— '^i  ~».  -  *■„  -\.  -  » 

>.V 


IT"  L’  P  ( N  +  I )  a  I  T  "  V  F  ( ’  } 
R  F  T  J  ’■’  \’ 

END 


'-,  ITT  VP.  C  1  ) 


SUBROUTINE  MARKO  () 

INTFCFR  J, K 
COMMON’  ./CALK/ 

CV'-'MCK  /TOFT/  TO (7 7,77) 

CO  ' ’  J  =  ’  ,  v 

CO  51  V  =  1 , X 

TO(J,C)=-1.0 
51  CONTINUE 

CV<J,J)=l.c 
57  CONTINUE 

CO  52  1 

TO (iron? (J)  , I  TOUR (J  +  l)  )=1. 7 
TC (I  TOFF (J+l) , I TOUR (J) )=1.7 
5  2  20^7  ls'r'^ 

TC ( ITCU? (!•*)  , I  TOUR ( 1 ) ) =1. 0 
TO  (I -CUP  (1  )  ,  I  TCP?  (V)  )  =1.  r 

P  f  T  VP  \‘ 

end 

S UB P " VT  17."  SCLCTY ( I , C, C,  N^R) 

I  NTFCF.P  I ,  Z,  ",  Y,  CO!.. 

1  NTL'CT  F  ? ,  3T\  (5  ;  ,  FTLA  (  5  )  , '  0(5)  ,»r; 


coy 

'  V  ~  \T 

/>.  r  LV 

v ,  i  to  up 

■*  < 
N_  •  * 

1  ''O  » ■ 

/TOFT/ 

10(77," 

r*  -'ll 
.* 

/ 1  N VT/ 

r  r  ( i  - j 

n  rss  ■ 

'  V  r  V 

/CIST/ 

rw  “>  /t  *7  ~f 
~  v  *  '  r  • 

c.  c 

7  =  1,  N 

7*0 

PC Y(J) =T 

C (I  TOUR 

TAP PAY  (7  ’) 


5”  ‘'  =  1,5 
F  =  ^ 

T  CM  p  =  ^Q,ir'>  n  n 

CC  51  7  =  3,"  V 

IF  (0  (IT  'VP  (2)  ,  I  TO  IT  (J)  )  .Cn.TEVF)  rCTCH 
I  F  ( T  '  F ?  A Y  ( T  TOUR  (  J  )  )  .CT.O.  ‘,)C'"'r^^l 
TFr  ?=C  (  I  T r •  ’ P  (2  )  ,  ITOUF  (.7)  ) 

F  =  J 

^  \'m  X  »•  .  ,  r 
•  i  i  .  „  I 

f  ^  /  7  v  ~ 


^  ‘T*  ,“v  <T  ^ 


•-•=”- 1 
o=5 


.  t  .  *  < 
2| 


V.  .»  .L.. 


gctcs  1 

G  4  TO FT AY  (ITOL’P  (  F)  )=l.f 

SY(K) =TFVP 

3  LA ( K )  *0  ( ITOL’P  (1)  ,  I  TCUP  (2)  ) -0  (TTOUP  (F-l)  ,  TT~UF.  (1)  ) 
ST'.'  (K)  =  I TCUP  (F) 

STLA (K ) = I TCUP  (F-l ) 

0  LA  (  K ) =3  LA ( K )  +1  ( ITll’O  (  F)  ,  I  TCUP  (F-l  )  )  -S  Y  ('<) 

n’KK)=f-i 

v  =y 

SO  tout  ii:ue 

STCP=n . 0 

IF  { .M.  EC.  ?)  <■  ”np=a(>oo  no 

T  EM  P  =  S  LA  (  l ) . 

CCL  =  1 

ic  s?  k  =  i , y. 

I F  (  S  LA(  K )  .  LF  .  TFM  P )  C  OTIS'*. 

TFMP=3  LA ( K ) 

CCL  =  K 

S3  CCSTIUUE 

C  =  3  LA  (  COL)  +STCP 
NCD=NO (COL) 

IT  (1*1  +  1) =STF ( CCL) 

IT (2*1+2)® S  TLA (COL) 

P  F.TUFS 
END 

SU3RCUTINF  ADJTFR  ( I  ,  P, SCO) 

1':tf.cfp  i,j,r,mco 
cc-‘vcn  /:■;?, l-</  s, iTOur  (71) 
c cry.?"-:  /toft/  Tn(7c,7.'>) 

C T-v  — ;  /T\VT/  IT(140) 

CO‘F’0*:  /TvTr /  I  'r  PTC0  (14  0) 

1C  *1  J  =  2 , N  C 1 

ITPTCF  ('100  +  2  -  J )  =  I  TCUP  (J) 

31  co':t  i;:ue 

10  P.2  J=2,’.10D 

ITCUF (J) = I T  PTCP (J) 

3  2  COST  IN’ IT 

T! (IT (2*1 -2)  , IT (2*1-1) ) =1 .  1 
TO  ( I T  (2*1-1)  ,  I T  ( 2  *  I  -2  )  )  *  1 . 

TO (IT (2*1-1) ,IT(2*I))=l.n 
TO ( IT  (2  *  I )  ,  I T ( 2  *  I  - 1 )  )=1. 2 
T 0(11(2*1)  , T T  ( 1 )  ) = 1 . ^ 

P.  =  3 

p  FTUP S 

rvn 


5  9 


1- 

K 

»* . 

i  ■  i 

m 


F 


i 


V'- 


0? 


91 


v.  • 


s  -7. ;;  -pc  i ':  e.  p  evt  r  r  ( ? ) 

I  N’T"  3FR  r 

~oyyc\  /.salk/  c,  iT-t1?  pi) 

ccvyo\’  /-'*  rr/  i t pt r °  (  i  i « . 

DC  ”2.  J  =  1,L’+1 

I  .  PTC?  (\'  +  ?-j)  =  i~c.r;F  (J) 
DDCTInuf 
DC  91  J=1 , M+1 

I  7CJ?  ( J  )  =  I  T PTC  ?  ( J  ) 
DCMTIVJF 
R  =  1 

RET'JP': 


Fl’NCTI 
REP  L  T 
C  ")v.  V.C": 
ix=i7i 
I Y= 1 72 
17*177 
IF( IX. 
I F ( IY. 
IF(  IZ  . 
TFV P=F 


CV  !'V7Cy(A) 

f;  -  f 

/??  ND/  I  v  f  I  y  ,  12 
VOCflX,  1  77 ) -  2  *  ( I  X/l  77 ) 

*VCP (IY, 17C)-2*(I Y/l 76) 

VCD  (IT,  1-7) -?*  (17/17=) 

LT .7)1 X*I X  +  2  7269 

LT.n) IY»IY+37277 

LT.  '»)  I7  =  r7  +  77  =  '>3 

LCAT  (IX)  /7  7 2  •'  9 .  ^  +F  L<~ -\ ?  ( I  Y )  /?  7  3  - 

*  *  V  C  C  (TFVr,  1 .0) 


7.  ■■  +FL-'NAT  (IT)  /=  7- 


R  FTC? y 
r 
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i  ■  ^AY,,ViV'..W 


rHiM.tr 


1/1  f 


b 


1  2 


1 

2 


0 

12 

11 


A£.gfHDIX  9 

LISTING  OF  HATRIX  UTILITIES  PROGRAM 


*  **  D  IS 

»-  * 

IMPLIC  I 

T  C 

i  niece  r 

I A 

COM:  MON 

/FI 

COMMON 

/NA 

WRITE  (* 

,1) 

WRITE  (* 

,2) 

WRITE  (* 

,3) 

W  R I T  E  (  * 

#4) 

WRITE  (* 

,5) 

WRITE  (* 

,6) 

W  R  I T  F  (  * 

,7) 

WRITE  (* 

,R) 

WRITE  (* 

,9) 

•  ,\ ' 


IX 


UTILITIES 
5  ( A ) 


0  (1  ',n,  122) 
V;AV  r  ( 1 ”2 ) 


WRITE (* , 12) 

READ (*,11)  12  MS 

IF  (  (IANS.  LT.  1)  .CP..  (  I  a  KS.GT.fi)  JCCT012 
IF ( IANS . EQ. 1 ) CALL  M*TPIX() 

IF (IANS. EC.21CALL  EXTPCT() 

I F ( IANS. EC. 3 ) C ALL  COPECT() 

I F ( I A  NS . EC. 4 ) CALL  MAT PET ( ) 

I F  (  IANS.  FQ .  5  )  CALL  r,IST() 

I F ( I A  NS  .  NF .  S )  C  CTO 1 2 
STCP 

FCRMAT(1X, 'This  is  a  set  of  utilities  to  be 
FORMAT (IX, ’  the*  orograx  ROUTE .  PUN .  They  ire 
FORMAT  ( IX,  '  create  and  manipulate  a  distance 

matrix* ) 


FORMAT (//IX, 
FORMAT (IX, ' 2 
FORMAT (IX,  '  3 
FORMAT ( IX  ,  '  4 
FOPMAT (IX,  '  ? 
FORMAT ( IX , ' 5 


used  wi th 1  ) 
to  he  Ip  '  ) 
"  a  t  r  i  x  . '  ) 


1 


•re. ate  a  o l stone 


1  matrix  from  a  laner  matrix') 


extract  a  sma 

correct  a  previously  store''  matrix') 
create  a  hardcopy  of  a  distance  matrix') 
list  names  with  their  row  number  in  matrix 
exit  this  program') 

F  ORMAT  (//l  X,  '  ".'.’hat  is  your  selection?  ’\) 

FORMAT (II) 

END 

SUBROUTINE  MATRIX () 

IMPLICIT  CHARACTER*  1  R  (A) 

REAL  01  ST 
I  NTI  CFE  N  ,  I  ,  J  ,  L 


r  VAN* 


'hi  ST/ 


C CM MON  /N AMR/  A NA M F ( 1 ^ 2 ) 
W  R I T  F  (  *  ,  3  D 
READ (*,32)  M 
WRITE (*,33) 


61 


,  v’  \  ‘  * 4  -** */  " 

*  “  .  •  *  .  ■  .  V  v“  v*  *  .*  <  ,  v*  \* 


CO  24  1  =  1,  N 

•'.'RITE  (*,35)  I 

PE AD (*,25)  A "AM  C ( I ) 

34  CONTINUE 

W  R  I T  F  (*,37)  m,n 
W  R I T  E  (  *  ,  3  ? ) 

DO  39  1=1, N-l 

D  (I  ,  I)  =^99. 99 

DC  49  J=I+1,M 

WRITE (*,41)  I,J 
READ (*,42)  DIET 
D (I , J) =D 1ST 
DO  IT?  L=1,1P9" 

IOC  CONT  I  t.'L'F 

D(J,  I)  =  D  I S  T 

IF (D 1ST. CE. 9) GOTO 40 

J=\’ 

I  =N  -1 

49  CONTINUE 

39  CONTINUE 

D (N , M) =9999 . 99 
CALL  STOP F (N ) 

RETURN 

31  FORMAT  (//IX,  '  How  larqe  is  the  oroblem?  ’ \) 

3  2  FORMAT (BN,  12) 

13  FORMAT (//IX, *  Input  names  for  home do  r  t  and  aids') 

35  FORMAT  (IX,  '  Number ',  13, '  Name?  '  \) 

26  FORMAT (A) 

3  7  FORMAT (IX,  13, '  RY  ',13,'  Matrix') 

3P  FORMAT ('  Input  distances  from  ooint  T  to  point  J') 

41  F OPM AT (  Distance  from  ot  ',13,'  to  nt  ',13,'  =  '  \) 

42  FORMAT ( BN , F7 . 2 ) 

END 

SUBROUTINE  FXTRCT ( ) 

IMPLICIT  CHARACTER*  1  5  (A) 

INTECEP  N,M,  I,J,  IFOWdEE) 

PEAL  S  D  ( 1  0  9 , 1  9 ) 

CHARACTER*  1  5  A  NM  (  1  r "' ) 

COMMON  /DIET/  D  ( 1  ^ R ,  1  "•'*) 

CC-'mcn  /:Dvr/  ?.  z  ( i  ^  ** ) 

WRITE  (*,1^) 

WRITF (* , 11) 

WRITE (*, 1 2) 

WRITE  (*,13) 

CALL  OBTAIN  (N) 

62 


.*<  VS  .*• 

•  -  -  .  ■' 

/.V.  -'I*.  ,'-  »-_V- W. 


VlsMs.s 

-  V  .  ■(»  m. 


< :  r  r  t  r 

u  r 1  c  ( * , :  - )  m 

DC  "?  ■'*  1  =  1 

a  R I T  E ( *,  2 1 )  I 

FORMAT  (IX,  n,  *  .  Row  number  fro-  cain  -ntrix:  ’\) 

R  FAD ( * , 15) IPCW( I ) 

CO  22  J  =  1  , N 

SD  ( I  ,  J)  =0  ( I  F.O/v  (  r )  ,  J) 

cont  INCF 

ASM  MM  (  I  )  =  A  N  A  M  F  (  I  RCV.'(  I  )  ) 

CONTINUE 
CD  2 J  =  1,N 

DC  2-1  1=1,  M 

D  (I  ,  J)  =SO(I  ,  I  °  O'. .  ( J  )  ) 

Z  O'-  N  T  T N  0  E 
0(J,J)=9on"'.r0 
A NAME ( J  )  =  A.  5V NM ( J ) 

CONTINUE 
CALL  .? TORE  ( '*) 

IT  URN 

FCFMAT (//l V, ’ To  create  a  scalier  matrix  fro-n  your  -nin1  ! 
FORMAT (ly,’  list  anco  -atrix,  enter  t'r  •  row  numbers  of) 
FnFv AT ( 17 , ’ "a  in  -at  t  ix  oils  which  yr  :  desire  in  the  ’ ) 
FORMAT (IX, ’scalier  ratr  ix*  ) 

FORMAT ( //! X, ’ ; h a t  is  the  site  of  the  smaller  natrix?  ’ \) 
FORMAT  (EM,  I  ;} ) 

!D 

'?  ROUT  IMF  0 'RECTO 
IMPLICIT  CMARA-'TFR*  1  5  (A) 

IMTFC-?  ",  T,  J,K,  I  '  *‘S  ,  IP,  IF 
P  L  '  L  •  I  c  T 

Zry"\Cy.  /~}l£7/  ~  ( I ’T1 ,  1  '  ^ ) 

C  OMMO  N  /N  A  VE/  A  NAM  F  ( 1  ”  ^  ) 

CALL  OPT' ON  (r: ) 

'..RITE  (*,14) 

/.'RITE  (*  ,  1  '  ) 

'..'RITE  (  *  ,  17) 

V.'PITF  {*  ,lf 
R  FAE  (  *  ,  1  a )  I  NS 
IF ( IANS, EC. 3) GOTO? 5 
I F ( I  * ’ ' S . F o . 2 ) OOT' 2 " 

IF ( IANS . NF . 1 1CCTC1 P 


M  ro 


sfa-(*,i:)  ic 

•■T-Trr.  (*  ,?o)  i  r ,  ic , d ( i p ,  ic) 

W  R I T  E  ( *  ,  2  4  ) 

R  FA 0  (  *  ,  1  ? )  DIET 
D  ( I  R  ,  I C)  =2  1ST 

do  is'1  :<=i,ior,r 

COOT  IMCF. 

0 (IC, I R ) =DI3T 
GOT  0 1 

WRITE (* ,2f  ) 

WRITE (* ,27) 

READ (*,12)  1= 

•W  R  I T  F  (  *  ,  2  3 ) 
p FAD (*,12)  IC 
DC  3  7  i  =  ir,::-i 
DC  31 

W?ITF(*,2-»)  I,J,C(1,J) 
WRITE (*,24) 

P  F  *  D ( * , 1 3 )  DIST 
IF  .DIST.CF.O.  i )  C-CTC  ^  2 
I  =N 


D ( I , J) =D  1ST 
DO  117  E=l,17a7 

COAT I  DUE 

D(J  ,  I)  =DIST 
CDMTIVCE 
I C* I +2 

D  (I  ,  I) =  ‘''>9?. 9 ° 

do::t  i  wit 

GOTO  1 S 

G-'.LL  3TCRF.  (I!) 

P  FTrjpv 

FORMAT (3M, 13) 

FORMAT (F7. 2) 

FORMAT  (//!’<, 1  This  utility  is  £: 
FORMAT  ( 1  X,  ’  i  "rtviouslv  stor:J 
F  ^  ?y  AT ( 1 M ,  *  1  -  Single  ‘-ntms 
EOF  V A  T  ( 1 X , ' 3  -  Exit  eo r  r  on  t  i  on: 
“rRMAT  (I  1) 


FOE  MAT (IX 
f op:-  a t  (l  y 
format  (IX 

F  ?!  "AT  ( 1  X 

Fr  I'M'  v 


' 7) w  n ur  be  t  to  c  o  r  r  o; 
'  Co  I  urn  n  number  to  no 
Tr?s.-nt  listener  fr- 
'"  :w  :  ist-mce  =  ’  V) 


:crr: - 
:a  f  i  1 


\) 

■t  ’M 


:•  t  i  n  -  ’  ) 


~u  on  t i 


*1  ~-nt 
n  :  *  '  ) 


tits’  ) 


' n t r y  of  -’-.M,  tyot 


:  --n  ->  1 1"- 


-**  fc.%  .*•  jO  _*•  m’; 


-  _  .  . 


F  3  p y a 7(!.v,  ' St irt ing  r o w  n u~> oe r  '  \ ) 
r  0 F  y  AT  ( 1 X ,  '  Starting  colon  narbrr 

■  v"i 


i  "B  POL’T  I  V Z  ATRRT  ( ) 

IMPLICIT  CHARACTER* 1“  (A) 

I NTFCER  I, J, M,  IS,  IF, JS , JF 
COMM  CM  ./CIST/  CirM’M 
COMMON  /MA*'F/  AVAME(1 
WRITE  (*  ,  ?.) 

'WRITE  {  *  ,  4) 

WRITE (* , 5) 

WRITE (*  ,6) 

CALL  OOTA.IV(N) 

CPF.M  (3  ,  FTLE  =  'VATRPT.  C?  TA  '  ,  STATUS  =  '\'  EK'  ,  FOPM=  ' ' 
15  =  1 
I  F=V 

IF (E.CT.  S'*)  I  F  =  5 ^ 

JS  =  1 
JE  =  9 

co  0  :=r?,i" 

IF(N.  LT.  I  "1)  J  F=K 
DO  15  J  = J S ,  JF 

WRITE  (3,  11)  0(1/0) 

CONTI MCE 
WPITF (3, 12) 

covt  i  A'  uf 

A  R I T  F (3/ 1 ?) 

js=jt+i 


J  F  =, 

JF 

+  0 

I  F  ( . 

*  V 

.  CT.M) 

Ji 

IF  ( 

T  c 

.  LF .  ‘.’) 

Cl 

II 

li. 

►-« 

I  E 

+  1 

I  E  = 

IE 

+  =  •> 

I  F  ( 

I  E 

.GT.V) 

I 

IF  ( 

IS 

. LF. M) 

r 

CLO 

SE 

(3) 

WRI 

TE 
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